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Frank A. Hagen, Doctor of Philosophy, 1976 

Thesis directed by: Frank B. McDonald 

Professor 

Physics and Astronomy 

The isotopic composition of cosmic ray Be, B, C and N has been 
studied using a new range versus total light technique. Special emphasis 
has been placed on the Be isotopes and, in particular, on the radioactive 
isotope ~®Be whose mean lifetime against decay (T^ = 2.2 x 10^ yr) makes 
it an ideal "clock" with which to measure the cosmic-ray age. The 
experiment consisted of a thin trigger scintillator, an acrylic pjastic 
Cerenkov detector and a spark chamber, followed by a totally active 
stack of 14 scintillation detectors. This stack of scintillators made 
possible the measurement of range, and also peirmitted the removal of 
interacting events by continuously monitering their identities along 
their trajectories through the. experiment . The experiment was c-..rri6'd 
by balloon to atmospheric depths ranging from 3.5 to 5.0 g cm"'^ residual 
atmosphere for a total e:<posure time of 23 hr- The flight was carried 
cut on August 15, 1973, from Thompson, Canada, 

Both the theoretical and empirical response of the experiment have 
been considered in detail. The data have been corrected for variations 
in detector response as a function of position and time, and also for 
variations in detector thicknesses. The mass resolution achieved by the 
experiment is given by ^ 0.047A. Comparisons are made with the results 
of other experiments after correcting the data for interactions in the 
detector, the varying energy windows (roughly 150 to 450 MeV/nuc) in 



which the experiment is sensitive to different isotopes, and production 

and destruction of nuclei in nuclear Interactions in the atmosphere 

above the experiment. In addition, the data have been corrected for the 

effects of solar modulation before being compared with model predictions. 

The data have been interpreted within the framework of the "leaky 

box" model of cosmic-ray propagation, where the cosmic rays are imagined 

to propagate freely within a box (e.g, the galactic disk) from which 

they have only a small probability of escape each time they hit a wall. 

The results indicate the survival of (55 +21) 7o of the in the 

arriving cosmic rays. In the leaky box model, this is interpreted in 

+6 b 

terms of a mean cosmic-ray confinement time given by = 5 x 10 

, + 1.0 _3 

yr, which corresponds to a mean density, n = 0.7 _q ^ atoms cm , of 

matter in the confinement volume. This result is apparently most con- 
sistent with the confinement of the cosmic rays to the galactic, disk. 
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CHAPTER I 


INTRODUCTION 

As understood at the present time, the cosmic radiation contains 
electrons, protons and heavier nuclei from throughout the periodic table, 
with energies ranging from the suprathermal to over a billion times what 
can be achieved using currently available accelerators. In the course of 
their nucleosynthesis, acceleration and propagation through space, these 
particles must certainly have participated in very high energy processes 
of various kinds, and therefore carry significant information character- 
izing these processes. The cosmic rays fill the galaxy to an energy 
density 1 eV cm"'^, which is comparable to that residing in visible 
light, the 3®K background radiation and the galactic magnetic fields. It 
follows that these particles must play an important role in the dynamic 
structure of the galaxy (Parker 1965, Badhwar et al. 1975). It is then 
of considerable interest to determine where the cosmic rays reside: 
whether they are confined predominantly to the galactic disk (Owens 1975, 
Dickinson et al. 1975), or perhaps fill a roughly spherical galactic 
"halo" surrounding the disk (Ginzburg et al. 1964, Suh 1974). The pos- 
sibility has even been suggested that the cosmic radiation may fill some 
extended region of extragalactic space (Setti et al. 1972, Brecher et al. 
1971, Sitte 1972). 

Of the various observable features of the cosmic radiation — compo- 
sition, energy spectra, temporal variations, and directional and positional 
anisotropies--the study of composition is perhaps the most promising for 
answering questions such as those posed above. For example, studies of 
the charge composition have revealed the presence of significant abund- 
ances of components (e.g. Li, Be and B) thought to be totally absent from 
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the cosmic -ray source. Such components must originate entirely in the 
nuclear spallation of heavier species. Thus the abundances of these 
species depend only on the amount of matter traversed by the primary 
cosmic-ray "beam". This leads to the conclusion that the cosmic rays 
have, on the average, traversed 

Ac = <rri>nirz^ ^ 5 

of material. Here <m> is the mean atomic mass of the material (gas) in 
the region where the cosmic rays propagate, n is the number density of 
atoms, V is their propagation velocity ( v ~ c is usually a good approx- 
imation) and is the mean propagation time corresponding to Xe* 

However, a knowledge of does not answer the question of where 
the cosmic rays reside (i.e. In the galactic disk, a galactic "halo", or 
extragalactic space). One parameter which is highly sensitive to the 
confinement region, however, is the number density of atoms, n: n ~ 1 for 
the galactic disk, n ~ 10 ^ for the galactic "halo", and n ~ 10“^ for Che 
interga lactic medium. Since X„ is known, however, the determination of 

C 

n and Tg are equivalent. As was first suggested nearly two decades ago 
by Hayakawa et al. (1958) and later emphasized by Peters (1963), this 
may be accomplished by a measurement of the abundance of the radioactive 
"clock" isotope . 

A direct experimental measurement of the abundance of ^^Be in the 
cosmic rays has long been delayed due both to the experimental problem 
of obtaining isotopic resolution and the rarity of ^*^Be relative to its 
neighbors. Attempts have been made in the meantime to approach the 
problem indirectly by the measurement of the charge ratio ^ (Shapiro 
et al, 1968, 1969; O'Dell et al. 1975). These attempts have suffered 
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from the rarity of Be with respect to the other isotopes of Be, in 
conjunction with the limited precision with which it is possible to 
calculate this ratio from theory in the cases of interest (i.e, the 
survival or decay of ^^Be). The first attempt at a direct measurement 
of the ^^Be abundance was by the University of New Hampshire group 
(Webber et al. 1973) using a balloon experiment floating at an atmos- 
pheric depth of ~ 3 g cm“^. This group has also recently reported 
results of a new, improved experiment (Preszler et al. 1975), Satellite 
observations taken outside the earth's atmosphere and magnetosphere have 
also recently been reported by the University of Chicago group (Garcia- 
Munoz et al. 1975, 1975a). The objective of the present thesis is to 
supplement and improve upon the above results by a measurement of the 
isotopic abundances of cosmic ray Be, B, C and N, with the primary 
objective of determining the cosmic-ray age by a measurement of the ratio 

lOfie 

Ee" 

Before discussing in detail the experiment and its results, it is 
desirable to sketch briefly our present understanding of the cosmic rays 
and the various theoretical models, particularly as these relate to 
studies of the composition of the cosmic radiation. The discussion is 
organized into three major parts following the temporal history of the 
particles: the nature of the cosmic-ray source (or nucleosynthesis’), 

the nature of the acceleration to cosmic-ray energies, and the nature of 
the propagation of the cosmic rays through space. With respect to each 
of these topics, we may inquire of the mechanism of the process, and its 
location in space and time. 

A. The Cosmic Ray Source 

Because of both the ultra high energies of individual cosmic-ray 
particles and the energy content of the cosmic radiation as a whole, the 
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cosmic-ray source has at one time or another been attributed to virtually 
every violent phenomenon known or postulated to occur in the universe. 

The earth, thunderstorms, the sun, magnetic white dwarfs, supernovae, 
pulsars, explosions of the galactic center, radio galaxies, quasars, 
and even the big bang itself have all been suggested at one time or 
another as cosmic-ray sources. Clearly, there is a strong need for 
experimental observations capable of restricting the options available 
to the theorists in constructing models. 

For example, the radio and optical emission of supernova remnants 
have been explained in terms of the synchrotron emission of electrons in 
magnetic fields (Ginzburg et al. 1964). For the particular case of the 
Crab supernova remnant, the emitting electrons are estimated to have 
energies as high as ~ 5 x 10^^ eV. Because of their overwhelming rate 
of energy loss by synchrotron emission, these electrons could not have 
been accelerated more than ~ 100 years ago. This is independent of their 
initial energy and far less than the age of the supernova, which was 
observed in the year 1054. Thus electrons are even now being accelerated 
to cosmic-ray energies in the Crab nebula. Accepting the rather uncertain 
hypothesis that the Crab nebula is a "typical" supernova remnant, and 
assuming that ~ 100 times the energy going into the acceleration of 
electrons goes into the acceleration of protons and heavier nuclei 
(which is the ratio observed locally in the cosmic rays), one may esti- 
mate the rate at which cosmic rays are being contributed to the galaxy 
by supernovae (Ginzburg et al, 1964), In particular, it is conceivable 
that a typical supernova dissipates ~ 10 erg in cosmic rays, which 
gives an average power of 3 x 10^® erg sec"^ if the frequency of super- 
novae is taken as one per hundred years. On the other hand, the ratio 
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of the volume in which cosmic rays are confined to their mean confine- 
ment time, — - ~ 10 cm'^sec"''', is characteristic of most cosmic-ray 

propagation models which confine the cosmic rays to the galaxy. For 
example, in the halo model of Ginzburg (1964), the halo volume is V 
2 X 10^® cm^, and the mean confinement time is ~ 3 x 10® yr. Assuming 
that the observed cosmic-ray energy density of ~ 1 eV cm~® is constant 
in time (which seems to have been the case for at least the last billion 
years; Van Loon 1973), we then estimate the power required of the cosmic- 
ray sources to be ~ 10^^ erg sec“^. This is in agreement with the above 
estimate of the power available from supernovae to within the accuracy 
of the numbers used. 

Such energy considerations have served to make supernovae and 
supernova remnants attractive sites for the acceleration of the predom- 
inant portion of the cosmic rays. This would not rule out contributions 
from other sources, and indeed there may well exist a hierarchy of dif- 
ferent types of sources both in our own galaxy and in others. Finally, 
we note that this result is dependent upon the existence of an efficient 
means of confining the cosmic rays to the galaxy, and thus retarding 
the rate of energy loss due to escaping particles. 

The characteristics of the cosmic-ray source can be probed by 
observations of the composition of the cosmic rays, as extrapolated to 
the source in terms of some convenient propagation model. For example, 
observations of the very heavy (Z ^ 30) component are apparently some- 
what suggestive of an r-process (as opposed to s-process) origin (Price 
et al. 1975). Since the r-process can only take place in an extremely 
hot environment with high density of neutrons present, this supports a 
model wherein the nucleosynthesis of the cosmic rays occurs in a 
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supernova explosion. The r-process does not (for example) occur as one 
of the normal cycles of nucleosynthesis in a stable star. Attempts have 
been made (with some degree of success) to calculate the observed cosmic- 
ray composition from models based on explosive nucleosynthesis (Hainebach 
et ai. 1975), which is a strong possibility as a mechanism for supernova 
explosions (Arnett 1969, Truran' et al. 1970). Such models have a number 
of free parameters, some of which may possibly be determined or consid- 
erably constrained by observations of the isotopic composition of the 
cosmic rays. For example, one such parameter is the neutron excess 
present in the source, which is strongly correlated to the abundances of 
and relative to ^^Fe. Another parameter, the temperature at 

which the nucleosynthesis occurs, would be reflected in the relative 
importance of various contributing nucleosynthesis chains, such as the 
explosive burning of C, 0 and Si, the CNO cycle, and so on. Each of 
these cycles produces certain characteristic isotopes, the observation 
of which would thus help to establish the temperature prevailing in the 
cosmic-ray source. 

B. Acceleration 

In the following section, the acceleration of particles to cosmic- 
ray energies is discussed. Features of the acceleration which may be 
reflected in the composition of the cosmic rays are of particular interest 
For these purposes, a detailed discussion of the physics behind the accel- 
eration is often not required, and therefore will not be given here. 
Features which' are of interest include the nature of the reservoir from 
which particles are accelerated, and how they are selected for acceler- 
ation. 

As suggested in the preceding section, supernova remnants are 
prime candidates for the cosmic-ray source. Furthermore, there is 
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evidence that the acceleration of cosmic rays is an ongoing process in 

supernova remnants (in particular, the Crab) . A potential means for 

accelerating these particles was provided by the discovery of pulsars 

(Hewish et al. 1968) — one of which happens to be located in the Crab 

nebula — and their interpretation as rapidly spinning neutron stars 

formed in supernova explosions (Gold 1968). The acceleration mechanism 

12 

would be provided by the enormous magnetic fields of ~ 10 gauss associ- 
ated with these objects. This is supported by the observation that the 
Crab pulsar is losing energy (as estimated from its angular deceleration.) 
at a rate essentially equal to that going into the acceleration of 
particles to relativistic energies (again assuming that ~ 100 times as 
much energy goes into protons and nuclei as electrons; Pacini 1973). 

It is implicit in all of the above that there may well be a sub- 
stantial time delay between the origin (or nucleosynthesis) of the cosmic 
rays and their acceleration to relativistic energies. Specifically, the 
nucleosynthesis would most likely take place in the supernova explosion 
Itself, while the acceleration appears to be a continuing process which 
might (for example) be directly attributed to the pulsar formed in the 
supernova explosion (e.g. see Gold 1968, Goldreich et al. 1969, Gunn et 
al. 1971), or might (on the other hand) be due to the action of rapidly 
moving magnetic clouds in the remnant (Ramadurai et al. 1972, Scott et 
al. 1975). Such acceleration would presumably occur after some degree 
of mixing of the debris from the explosion with the interstellar medium, 
which would be reflected in the cosmic-ray composition (Reeves 1973) . 

The existence of a time delay between the nucleosynthesis and acceler- 
ation of the cosmic rays can be tested by measurements of the abundances 
of certain isotopes (e.g. and ‘^'^Ti) which decay exclusively by 
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electron capture, and are presumed to be produced by the cosmic-ray 
source (Shapiro et al. 1975, Soutoul et al. 1975). At rest, these 
isotopes are neutral and will decay at characteristic rates. Once 
accelerated, however, they are stripped of all their atomic electrons 
and thus become stable. We thus have radioactive clocks which are 
turned on at the time of nucleosynthesis and conveniently turned off 
when the actual acceleration occurs. Although isotopes which decay 
exclusively by electron capture are best for these purposes, it is pos- 
sible that one may also use isotopes' with alternate modes of decay (e.g. 

56 59 

^Ni and Ni), provided the alternate lifetime is sufficiently long 
(i.e. > a few million years, depending upon one's model for propagation) 
to ensure their survival in the cosmic rays. 

Models have also been proposed in which the acceleration occurs 
virtually at the same time as the nucleosynthesis, i.e. by relativistic 
shockwaves generated by the supernova explosion (Colgate 1969). How- 
ever, such models suffer from a number of difficulties. One questions, 
for example, the survival of the heavy nuclei in the extreme environment 
of a relativistic shock wave. It also appears inevitable that the intense 
streaming of the cosmic rays away from their source when created in such 
intense bursts will set up self-generated magnetohydrod^mamic waves 
(Wenzel 1973). These waves then serve to efficiently absorb energy -by 
scattering the cosmic rays. Such a result implies that the cosmic rays 
would escape with at most a few percent of their original energy. This 
drastically increases the power requirements on the cosmic-ray sources, 
probably to unacceptable levels. Also, the dissipated energy would heat 
the interstellar medium, again to a degree which may well be unaccept- 


able (Wenzel 1973). 
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The remarkable resemblence of cosmic-ray abundances to those 
characterizing the solar system (which may also be explained, in large 
part, by models based on explosive nucleosynthesis) has been noted 
essentially from the first observation of the nuclear component (Bradt 
et al. 1950) « The correlation is considerably enhanced by taking account 
of ionization cross sections (Casse et al. 1975), which suggests that the 
cosmic rays may result from the selective acceleration of material char- 
acterized by essentially solar system abundances. The selection mechanism 
would be ionization--i-e . only ionized atoms are accelerated. Such a 
model would, of course, predict no anamalous enhancements or deficiencies 
of specific isotopes as compared to other isotopes of a given element, 
as the ionization cross sections are alL the same. This line of reason- 
ing can thus be easily tested by measurements of the isotopic composition 
of the cosmic rays. 

C. Propagation 

Having discussed the origin and acceleration of the cosmic rays, 
one still needs to consider their propagation in interstellar space 
before being able to make comparisons with experimental observations. 

Most current models for the propagation of cosmic rays are based on the 
diffusion model originally proposed by Ginzburg et al. (1964). Although 
cosmic rays are constrained to closely follow magnetic field lines in 
this model, they may be scattered back and forth along the lines by 
irregularities in the lines. Furthermore, the actual orientation of the 
field lines in space is considered sufficiently random that a small 
drift (or scattering) of a particle across the field lines effectively 
randomizes Ics trajectory through space. Thus the propagation has been 
approximated as a diffusive process, with the actual diffusion coefficient 
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used left as a free parameter to be adjusted (within limits) for the 
optimum agreement with observations. 

If the characteristic pathlength required for a particle to 
escape in a diffusion model is much greater than the typical dimensions 
of the confinement region, then the diffusion model may be quite well 
approximated by a "leaky box" model (Jones 1970, Owens 1975a). In this 
model, the boundaries of the confinement region are pictured as the 
sides of a "leaky" box, which have a high efficiency for reflecting 
particles back into the box. The fundamental parameter in such a model 
is the probability of leaking out of the box. This is given by the 
leakage pathlength, which is related to the mean time for escape 
from the box, Tg, by equation I.l. If there are radioactive isotopes of 
mean decay time T^j ~ Tg involved in the problem, then the model also 
becomes directly sensitive to the mean escape time, Tg, and thus to the 
number density, n, of atoms in the box. The parameter can be deter- 
mined if we assume that certain components of the locally observed cosmic 
rays are absent in the source. Their origin must then be entirely in 
the fragmentation of heavier species in the interstellar medium, a process 
which can easily be modeled. For example, Li, Be and B are all easily 
destroyed in hot stellar interiors, and for this reason are thought to 
be absent in the cosmic -^ray sources. A calculation based on these 

elements gives 5 g cm ^ (Shapiro et al. 1973). Assuming a density 

-3 

n ~ 1 hydrogen atom cm in the confinement volume (which would be con- 
sistent with confinement of the cosmic rays to the galactic disk) we 

£ 

then obtain Tg =>* 3 x lO'^ yr. from equation I.l. Alternatively, we may 
adopt a "halo" model in which the cosmic rays are confined to a quasi- 
spherical volume surrounding the galactic disk. In this case the gas 
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dens4.ty (n 2 : is much lower than for the disk, so that a much 

longer escape time (i.e. Tg 2 ; 3 x 10® yr) is required for consistency 
with the observed leakage pathlength. It is thus apparent that an experi- 
mental determination of the mean escape time, Tg, would be a sensitive 
test for discriminating between these models. 

As has already been pointed out, an experimental determination of 

the mean escape time, Tg, is possible by a measurement of the .abundance 

of a radioactive isotope which is completely secondary in origin (i.e. 

absent from the cosmic-ray source), and also has an appropriately long 

mean lifetime against decay (i.e. ^ Tg). Several such isotopes have 

been suggested, including (T^j = 2,2 x 10^ yt), (Tj = 1.2 x 10^ 

yr) and (T^ -Ax 10^ yr) ; in extreme cases even the isotope 

9 

(T (3 = 1.9 X 10 yr) may be of use. Most attention so far has been 
focused on ^®Be because of its clear secondary origin, its nearly ideal 
lifetime (for a disk confinement model), and also because it is experi- 
mentally the most accessible. The experimental problem of separating 

(especially from the more abundant neighboring isotope ^Be) is still 
very difficult, however. For this reason, most of the effort so far has 
been concentrated on the indirect approach of measuring the charge ratio 

.5®. (Shapiro et al. 1968, 1969; O’Dell et al. 1975). Since decays 

B 

to this charge ratio is the most sensitive of any to the question 

of the survival of ^*^Be in the cosmic rays. The precision obtainable by 
this approach suffers greatly from the very low abundance of ^^Be relative 
to the other Be isotopes, in conjunction with the limited precision of 
the theoretical predictions which can be made. Nevertheless, the 
analysis has been performed, with results which are most consistent with 
the' complete survival of ^^Be in cosmic rays of energy E ^ 1 GeV/nuc 
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(Shapiro et al, 1968, 1969; O'Dell et al. 1975). The first attempt to 
directly measure the abundance of ^^Be in the cosmic rays was by Webber 
et al. (1973) , who used a balloon borne experiment floating beneath a 
residual layer of ~ 3 .g cm"^ of atmosphere. This group has also recently 
reported results on a new, improved experiment (Preszler et al. 1975). 

The only other published results are those of Garcia-Munoz et al. (1975, 
1975a), who report results from the IMP 7 and IMP 8 satellites. While 
Garcia-Munoz et al. feel their results are consistent with the complete 
decay of ^^Be (i.e. a very long leakage lifetime), the results of 
Webber et al. appear to indicate the survival of a significant fraction, 
of the ^^Be. 

It is the objective of the current thesis to supplement and extend 
the above results by a measurement of the isotopic abundances of the 
elements Be, B, C and N, with special emphasis on the question of the 
degree to which ^^Be survives in the cosmic rays. Our results indicate 
the partial survival of ^^Be, and result in the estimate Tg = 5^ x 
10^ yr for the mean escape time of cosmic rays. This is entirely con- 
sistent with models which confine the cosmic rays to the galactic disk. 

It is not consistent with models in which the cosmic rays are confined 
to a spherical halo region, but may be consistent with a hybrid model in 
which the cosmic rays are confined to a considerably flattened halo region 
surrounding the galactic disk. 
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TECHNIQL'KS FOR THE ^XASUR:-:M^^•T OF 
A. Tho Devclopr.cnt of Expor i=vnta 1 Tt*chniqucs 

Boloro discussing in detail the general problcr. of determining 
the mass of a cosnic-ray particle (section B) , and Cl;e specific expi?ri- 
mental approach of this thesis (section C), a brief description yill 
first be given of the deve lopmt'nt of our kncvledgc of the cor.pcsition 
(chemical and Isotopic) of the cosmic radiation. Wc begin with the dis- 
covery of the existence of a significant component consisting of nuclei 
heavier than protons (Frier et al. 1948), which followed closely upon 
the association of the predominant component with protons by its highly 
penetrating nature (Schein ot al, 1941) and positive ctiarge (Johnson 
1933, Alvaros et al. 1933, Rossi 1934). Although the existence of a 
heavy component was soon cor.firtned by a number of other experim.cnts (e.g. 
see Singer 19S8 and references therein), it was to be a decade before 
experimental resolution was improved to the point where a consensus began 
to emerge or. the detailed charge cor.po«;it ion of the cosmic rays. In 
particular, there were s ign i f ic.ant discrepancies in the measured .ibund- 
anecs of the elemi»nts Li, Be and B as well as the odd 7, nuclei (e.g. S, 

F, Na , Al and P) , which tend to have abundances considerably reduced from 
tho.se of their even 7. neighbor:*. At the present time, exper iavnta I tech- 
niques h.nve improved considerably, to the point where excellent charge 
resolution has been obtained by groups working at several laboratories 
for elements with 7 t; 26. 

1. Basic Detectors . From the beginning, there has been a wide 
variety of detection schemes applied to the study of the cosmte-rny 
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composition. These have included both passive detectors (e.R. era.jlsions 
and plastics) and active detectors (c.g. cloud chanbers» spark cl)ari>ors, 
proportional counters, scintillators and solid state detectors). Although 
the passive detectors dominated at first, rapid Improvements in the 
resolution and areas of the various active detectors enabled them to 


overtake the passive detectors. 

2. . Farly investigations of the chemical conpos it ' ."sn made 

use of the fact that the georiagneCic field excludes low energy particle;: 

from approaching the earth at low latitudes. 3y going to su f f i c icnr 1 y 

low latitudes, one was therefore able to select highly re l.it ivlstic 

particles (F. > 3.5 CeV/nuc) whose ionization energy loss is propnrtion.il 
2 

to Z . In thi.s way. It was possible to determine the cliarge compn;: it ion 

of the cosmic radiation at high energies by a simple measurement of Che 

dE 

ionization energy loss or stopping pewer, 

3. ^ X C. An imporcar.r step towards extending r^-asurer-ynt s 

of the chemical cempos it ion to particles of arbitrary energy wa;: taken 

by Webber and McDoivild (1955), who added a measurement of Cerenkov radi- 
d“ d^ 

aClon to the ^ ne.isurement . Since both ^ and the em.l::sicn of Cerenkov 

radiation depeml only on 7. and B, this .illywed the determinati<;n of both 

these paranc; ters by the s itni 1 taneous solution of two equations (i.e. one 

for ^ and one for the Cerenkov radiation). In tills way, the chemical 

composition of the cosmic rays may be mie;asured over a wide r.mge of 

2 5 

energies. Tliis range of energies (^- 10 to 10 HeV/nue is essentially 
defined by the threshold and saturation veloc;ities characteristic of the 
available radiator materials. 


4. •— X F. The slnxil ta neons mc.asureoHnt of the ion i /.at ic:: energy 

dx 

loss and the total kinetic energy .of a particle provided a s«-cotul 
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technique CAp.ibIc of seasuring the chec:lcjl cocposlClon as a function 

of energy (U\-bbcr et al. 1972). The basic principle of this approach 

Is slcilar to that of the technique described above, with the Cerenkov 

ocasurer.enC replaced by a total energy ceasurcc:cnt , The difference is 

that the energy depends on taass (A) and velocity (P), whereas the 

Cerenkov eoission of a particle depends on Z and g. One is therefore 

required to assure a one to one correlation between charge and mass 

(c.g. A =: 2Z for Z 20) in order to complete the set of equations for 

the unknwnn A, 7. and g. In one form or another, this technique has 

been applied over a range of energies extending frors ~ 10 keV/nuc up to 

~ 100 CoV/nuc. The liir.itation at lower energies derives froa the 

dE 

racroscopic thickness of tf-.e ^ detector, while that at high energies 
results from problems associated with the construction of a practical 
detector capable of cornpleteLy absorbing a particle's energy. As one 
goes to energies E > 100 MeV/ruc, it bcc;xr.es increasingly Ir.probable 
that a given particle will stop in a detector without undergoing .i 
destnictive nuclear interaction. Although it in still possible to 
estiryite the energy of a p.irtlclc under sue!; circumstances, the resolu- 
tion of the r.ea si:rerer,t is considerably reduced. 

dE 

As resolutions of E and ^ detectors continued to improve, even 
Isotopic resolution becar.e possible by this lech.niquc. Tr. this case, 
the assumption of a one to one correspondence between A and 'L is replaced 
by the more realistic assur.pt ion of discrete valued charge and mass dis- 
tributions (see .section 3 for a detailed discussion of the analysis). 

The mass resolution obtained by this technique characteristically 
deteriorates as one goes to heavier particles. tie and Me vrere thus 
the first isotopes separated. At the present tine, the isotopes resolved 
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by this technique extend lo nitrogen O^ebber et ol. 1973, 1973 j) ,ind 
beyond (Creimir 1972, Mewaldc <*i al. 1975). 

5. Magnetic Spec trore ter s . Recently, large superconducting 
cagnetlc .spec irorcc ter s have al^o been applied to the study of the charge 
and Isotopic coeposir. ion of the cosrjic rays (IXauber 197:). tliesc 
Instrururnt.s alloc? tt;e detemlr.at ion of both, the sign and aagnitude of a 
particle’s ch.irge, and thus have been u.sed to search for antisuittor in 
the cosn^ic rays (Bu f f ingt.on et al. 1973, Golden ct al. !97^»). It has 
also been possible to use the geomagnetic field as a rigidity filter 
(Jullusson 1975 , Dwyer ot al. 1975) in order to study the i.sotopic 
composition of the cosmic rays at energies of a few CeV/nuc , which are 
very much above the energies accessible by other -.ethods (i.e. K 'h 

500 MeV/nuc ) . 

6, Tra lec tory De f in it Ion . As geometrical factor:; (i.e. active 
areas) continued to increase in attempts to measure ever rarer particles, 
some means for trajectory definition of particles in oxpcritw'nls became 
essential. Ttii.s has been accomplished using hodoscoper, (thin strips of 
scintillator or flash tubes, for example), optical and digital spark 
chambers, and multiwire proportional counters. Using such detector:;, it 
became possible to correct for the zenith angles at which particles were 
incident on detectors, and also to correct for variations in detector 
efficiency over their surfaces. 

The various possible experimcnt.il approacl'.es to the determij-.at ion 
of the isotopic composition of the cosmic radiation will be di::cusr,eC in 
greater detail in s-ctlon B, along with the approach actually adopted for 
this experiment (section C). 
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B. The MeasuritcxTC of >lass 

For our purpose:;, an in<!ivi<iual cossic-r^y p^irclcle rany be 
described by throe parnneter:; ; its nwi:»s (A), charge (7,) and velocity 
(S). Tl^is description ignore:; properties such as sense of v«»locily, 
degree of ionization and levels of nuclear and atomic excitation, but 
will suffice for the pre:ient discussion. As it turn?; out, the charge 
may sor.etirn'S he determined with no a priori knowledge of the other two, 
while velocity may be determined knwlng only the charge. The deter- 
mination of the rjjss, on the other hand, re<;uires a knowledge of both 
charge and velocity and, in this sense, is the most difficult r^-asure- 
oent of the three. One has available essentially five observable:; from 
which to determine the mass; the range of the particle in some absorber, 
its energy, its rate of energy loss to ionization, its rigidity and its 
emission of Cerenkov radiation upon passing through an appropriate 
radiator. Since there are only three parameters to be determined (A, 2 
and B), the problen is over specified, and it is sufficient to treasure 
a subset consisting of only three of the five observables (assuming, of 
course, that they are all Independent). In fact, it is soc.otirars suf- 
ficient to measure only two of the five by aiwuiming the particles one 
will observe are stable against spontanecnis nuclear decay, and thus have 
A 22 . In the following five subsections, each of the five observables 
will be discussed in turn. 

1. Energy . The kinetic energy of a particle of mas:: A, charge 2 
and velocity % is given by 

= (r-ij = A\siA) (n-u 

where y denotes the relativistic Lorenz factor. As will be customary 
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for Che rccalndcr of this thesis » nass number and charge are Indicated 
by superscripts and subscript.'!, respectively, on the left hand side of a 
variable. The first equation ru^rely states the usual relativistic 
fonnula for the kinetic energy, vhilo the last equality give.s the scaling 
In A and Z. The dependence of on 0 Is shown in flgtire II. 1. 

The energy of a particle i.s often measured using a ca lor It&et r i c 
approach in which the particle’s energy l.s used to ionize atre-.s in soov' 
suitable absorber. Tlic energy deposited Is then determined by sor/Y means 
characteristic of the specific type of detector being used. ror example, 
the ionization calorimeter tzay consist of scintillation counters, in 
which case one observes the light emitted in the rccorbinat ion of 
electrons and ions. Alternatively, solid state detectors may be used, in 
which case one observes .in increase in the conductivity of Che dctectcjr. 
Finally, one may use passive dctector.s such as stacks of ertilsions or 
thin plastic sheets, in which ca.se one looks for actual dansige done to 
the material. In some c.ises, passive absorbers are inserted between 
layers of active rvuerlal (c.g. scintillators), but this is done at tb.e 
expense of the precision of th<- measvirem(.*nt . The resolution of such j 
mcasurenent l.s physically limited by tlu> Poissonian ;; t.i 1 1st ic s assoc i.ned 
with processes such as Che production of ion-electron pairs, tiu’ con- 
version of such p.iirs to photons and the conversion of photor.s to photo- 
electrons in 3 photomultiplier tube. For example, If Wp^^ i.s the energy 
deposited per photoo lectron produced ,. the resolution is giver, by 



For scintillation detectors, 20 keV/photoe Led ron so 

resolution better th.in 17. may be obtained for K > 200 McV. 


(II. 2) 

that energy 
For other 
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PROTON ENERGY AND 
GRADIENT VERSUS VELOCITY 



Figure II. I . The energy of a proton and il*; velocity 
derivative as functions cf velocity. 


E (MeV) 
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types of detector (such ns solid state), nay be considerably smlier, 

and the energy resolution is accordingly increased. Practically, !;uch 
resolutions are often far in excess of the requirements and siot fully 


realized in laany experiments. 

2. Rate o f EnerRV l.oss to lonizat ion . The rate at which a 
particle loses energy by ionization, or stopping power, is given (in 
the Born approximation) by the Bethe-Blocli formula (Bohr 191!), Beth.e 


1930, Bloch 1933) 


Vil) = 'V Zr n^ff/ 

^ (11.3: 

^ 'm 

This expression breaks d<rwn .at very high energies, where t!.e In y 
Increase of rhe term, in brackets is eventually checked by the ’'density 
effect" (Sternhe Ir.er 1952, 1961). ror particles of very low velocity, 
on the other hand, electron pickup becomes important and 7. rxist be 
replaced by the effective (tine averaged) charge of the particle. 

Neither of these effects are important for our considerations. Again, 
the scaling of the Ionization energy loss In mass and charge is given by 
the last equality of etjualion 1 1. 3. The dependence of the proton ioni- 
zation energy loss, velocity, g, is shown in figtire I 1.2. 

A particle's rate of energy loss by ionization may be measured 
using proportional counters, ionization counters, thin scintillation 
counters or solid state detectors, gnralsions and chin sheets of plastic 
may also be used. Resolutions of the order of ~ 17, (dominated by Landau 
fluctuations; Symon 19^8, Rossi 1952, Tschalar 1968, 1968a) are obtaln- 
oblc, depending upon the thickness of the detector, and the nature of 
the primary particle. 


(11.3) 
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3. Range » Range Ik related to the stopping power, and energy, 
2 

E “ ( 7 -l)AEipC , according to 


W)- 


. /I 


^ I J£ 




= m 




where the scaling law given by the last equality follo-^s from that for 

dE 

the stopping power, rigure il.3 shows the velocity dependence of 

the proton range, Jr. 

The range of a particle is best measured by noting where the 
particle stops in a stack of emulsions or thin plastic detectors. More 
crude measurements using (thicker) solid stale detectors or (even thicker) 
scintillation detectors have also been of some use in prt'vious'. experi- 
ments. The range of a particle fluctuates due to t!ie stochastic nature 
of the energy loss process. Tlie distribution is characterized jpproxi- 
mately by a Gaussian with standard deviation 0.0lRj\ ‘ (^ternheimer 

19f>0), and imposes a physical licit on ti;e precision with which a rx>an- 
Ingful measurcovent can be mde. 

Rigidity . Rigidity, the fotirth observable, is given by 




(11.5) 


Once again, the last equality gives tb.c scaling. The dependence of the 
proton rigidity, [p, on velocity is sl.own in figure 11.4. 

The measu recent of the rigidity of a particle is accomplished by 
measuring its radius of curvature in a known magnetic field. This makes 
use of the relation 


p 3_ILL 

2e c^j 6 


( 11 . 6 ) 


where $ represents the angle between the particle's trajectory an<i ib-e 
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oagnvtic field. The gyror.idius , , i:s determined using position 

sensitive detectors, such as cloud ch.imbers, bubble ch.imbers, rulti- 
wlre proportional counters, spjrk chambers, emulsions or plastic detector 
Spark chambers and oailtiwire proportional counters have the advantage of 
being corpatible with digital readout^ while the others (also optical 
spark chambers) are generally capable of higher spatial resolution. 

5. Cerenkov R adiation , finally the Cerenkov response (in units 
of the response to a p ■ 1 proton) is given by (Jelley 1958) 




(II. 7) 


n»o last c(;ualiLy again gives the scaling, and the proton Cerenkov 
response, Jc, is illustrated in figure 11.5. 

Cerenkov radiation n-ry be observed by placing a radiator wit!, a 
suitable index of refraction n (depending upon the range of velocities 
one intends to study) in a diffuse reflecting box. The liglil emitted 
by the radiator is then observed by photoru ! t ip 1 ie r tubes, Titc dor.iiwr.t 
contribution to the resolution obtained is the statistical process of 
collecting photons and converting them to photoc lec t rons , w!;ich is 
governed by Polssonian statistics. If a Cerenkov response of unity cor- 
responds to € phot oc lect rons , then the Cerenkov resolution is giver, by 

on = (^) (II.8) 

Values € 10 are possible, depending on experimental details such as 

the thickness of the radiator, efficiency of plioton collection, and so on 
6. M. 1 SS Determination . Given observations of some subset of the 
five observables discussed above (i.e. E, S, P and C) , the relevant 
subset of equations II. 1, II. 3, II. 4, II. 5 and 11.7 may be solved 
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simultaneously for the mass of a particle, together with its charge and 
velocity, We now consider the question of what is the optimal subset of 
the observables which one can observe. Suppose, for example. Chat we 
have a ^%e incident on an apparatus with velocity p = 0.5. We have 
shown in figure II. 6 how the analysis might proceed if we measure (for 
example) the observables E, ^ and/or R. The appropriate numerical 
(i.e. observed) values have been substituted into each of equations II, 1, 
II. 3 and II. 4, and the resulting relations solved for the mass. A, as a 
function of charge and velocity. The results are illustrated in figure 
II. 6 for the particular cases where Z = 3, 4 and 5. Consider first a 
two parameter analysis utilizing the observables E and R. The various 
possible results of such an analysis are represented in figure II. 6 by 
the intersections of E and R curves of the same charge, Z. In particular, 
the solutions are 




A 


M) 


) 



(II. 9) 


for the present case. It is apparent from the above example that a two 

parameter analysis sometimes may be sufficient to uniquely determine all 

three of the unknowns A, Z and P when supplemented with the requirement 

that the observed particle be one of the known, stable isotopes. In 

particular, for the present case this immediately rules out the solutions 
20 . 5 

Li and B. The results of a three parameter analysis may, in similar 

fashion, be represented as the simultaneous Intersection of all three of 

the appropriate curves on a diagram such as figure II. 6. In table II. 1, 

the results of the various possible two parameter analyses involving the 
dE 

observables E, R and P have been tabulated for Z = 3, 4 and 5. We 
have again taken the actual event to be a ^^Be of velocity p = 0.5. 
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Figure II. 6. Mass as a function of assumed velocity derived 
from measurement of energy, ionization energy 
loss or range of a ^^Be with velocity g = 0,5. 
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TABLE II. 1 



Results 

of Various Possible 

Two Parameter Analyses 


Involving E, 

^ and R for Be 

Event of Velocity p 

= 0.5 

Ana lysis 





Mode 

Z 

A 

Sl 

Stability 


3 

21,20 

0.360 



4 

10.00 

0.500 

Stable 


5 

4.85 

0.650 


ExR 

3 

20.00 

0.372 



4 

10.00 

0.500 

Stable 


3 

5.00 

0.642 


ExP 

3 

5.30 

0.630 

Marginal 


4 

10.00 

0.500 

Stable 


5 

16.45 

0,403 



3 

22.40 

0.360 


4 

10.00 

0.500 

Stable 


5 

4.80 

0.650 


^xP 

dx^ 

3 

11.10 

0.360 


4 

10.00 

0.500 

Stable 


5 

8.40 

0.650 

Marginal 

RxP 

3 

8.50 

0.464 

Marginal 


4 

10.00 

0.500 

Stable 


5 

11.25 

0.540 

Stable 
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Since the event considered would be below Cerenkov threshold, only an 
upper limit on the velocity (i.e. p ^ 0.67) can be extracted from a 
Cerenkov measurement. For this reason, the observable C has not been 
considered in table II, 1. From the table, we see that the only reason- 
able result (i.e, one yielding a stable particle) of a two parameter 
analysis involving combinations of the observables E, and R is a 
of velocity p = 0.5, as desired. Furthermore, we gain little more than 
redundancy (which may be useful in the identification and elimination of 
background) by combining all three of these observables in a three para- 
meter analysis. In particular, the degeneracies present in the various 
possible two parameter analyses (without the stability condition) are not 
removed by the three parameter analysis unless the measurements are 

extremely precise. This is readily apparent, both from table II. 1 and 

dE 

figure 11,6. The two parameter analyses using either E or in com- 
bination with P result in extraneous solutions which are only marginally 

rejected (or accepted) by the stability requirement. A three parameter 

dE 

analysis using all three of these observables (E, ^ and P), however, 

easily removes the degeneracies (even without the stability condition), 

dE 

as opposed to the case where E, ^ and R are used. Finally, the two 

parameter analysis using R and ? is the most ambiguous. None of the 

extraneous solutions with 2^Z^8 generated by this approach can easily be 

rejected by the stability condition. In this case, a third parameter 
dE. 

(either E or is absolutely essential for the analysis. We may con- 

dE 

elude from the above discussion that any combination of E, and R is 
adequate for a two parameter analysis. On the other hand, the combination 
E, and P yields the best three parameter analysis. We also note in 
passing that both of the observables E and R require that a particle be 
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stopped in some absorber for their determination. The probability of 

being able- to do this without having the particle destroyed in a nuclear 

interaction decreases exponentially as the range of the particle (and 

also its energy) is increased. As a result, an analysis involving the 
dE 

measurement of ^ and P (which do not require that a particle be stopped) 
becomes increasingly attractive. The measurement of Cerenkov response, 

C, is similarly attractive for velocities above the Cerenkov threshold 
(i.e. p 2 : 0.67 for a radiator with n = 1.5). 

7. Resolution . The above discussion of various analysis schemes 
for mass determination is essentially independent of the resolutions 
obtainable or other technological constraints. We now consider the 
resolution with which it is desirable to measure the various observable 
parameters in order to achieve isotopic separation. Consider, for example, 
a mass distribution of the form 

This represents the result of measuring two neighboring isotopes (of 
masses = 0 and A 2 = 1) with abundances in the ratio k:l, respectively, 
with a mass resolution given by cr ■ For the purposes of this discussion, 
the mass scale is arbitrary to within a constant, i.e. the isotopes 
could as well be of masses = 9 and A 2 = 10, in which case one need 
only make the substitution A -♦ A-9. The condition for the existence of 
an extremum (peak or valley) in this distribution at A = A„ is 



( 11 . 11 ) 


If we choose k = 2 (which is approximately what we expect for the case 
9 10 

of cosmic ray Be and Be), the extremum location, A^ , is related to 
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the resolution, as shown in figure II. 7. In this figure, the value 
of the distribution function at the extremum (i.e. the height of the 
peak or valley) is also indicated. For example, consider the case where 
o = 0.5. From figure II. 7, we see that there is only one extremum (at 
= 0.085) in this case, which must correspond to the major (A^^ = 0) 
peak. The minor (A 2 = 1) peak is apparently lost in the tail of its 
neighbor. If we now Improve the resolution to O' = 0,3, there are three 
solutions to equation II, 11. These are at = 0.005 (the major peak), 

= 0.599 (the valley) and A^ = 0.991 (the minor peak). We also see 
that the height of the distribution is f(0,005) = 2, f(0,559) = 0.674 and 
f (0,991) = 1,10 at the major peak, valley and minor peak, respectively. 

The condition on the mass resolution for the appearance of a 
valley between the two distributions is cr ^ 0.35 in the particular case 
where k = 2, as may readily be seen from figure II. 7. This result may be 
generalized to arbitrary k by writing s '^crit* '^crit ^^epresents 

a maximum of the function O' (a^) defined by equation II. 11. Again, this 
is readily apparent from figure II. 7. For the general case, we thus 
require 
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( 11 . 12 ) 


' i-kJ 

where the second equality is merely a restatement of equation II. 11. 
From equation 11.12, we' see that a knowledge of any one of the three 
parameters k, and defines the other two. In particular, the 

dependence of shown in figure II. 8. This figure shows 


how the requirement on the resolution for the appearance of a valley 
becomes more strict as the relative abundance of the minor isotope 



t(Ao) AND <r(Ao) 



Figure II. 7. Resolution as a function of extrema locations for a bi-Gaussian 
distribution with amplitude ratio k = 2. Values of the distri- 
bution function at the extrema are also indicated. 
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Figure II. 8, Critical resolution required for the 
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distribution as a function of amplitude 
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(e.g. Be) is decreased. 

We emphasize that the existence of a valley between two neighbor- 
ing isotopes in a mass distribution is only one of many possible defin- 
itions of the critical resolution required to "separate" the two isotopes. 
One may even rely completely on deconvolution techniques in which case the 
"separation" criterion Is on the accuracy to which one knows the resolution 
and indeed the shape of the mass distribution, rather than on the magni- 
tude of the resolution itself. 

In order to estimate the resolution obtainable by the various 
analyses described above, we formalize equations IT.l, II. 3, II. 4, II. 5 
and II. 7 as follows: 

- Oi <11. 13) 

A 

Here 0. represents the observed value of a given observable, while zO,* (8) 
gives its functional dependence on Z, A and p. For example for the 
observable E, we have 




and 


(11.14) 


(11.15) 


Equation III. 13 may be solved for the mass, A^(p), as a function of 0^^^, 

Z and p. The functions Ag(p), A^g(R) and Aj^(6) are shown in figure II. 6 

dx 

for the particular case where the observed particle is a Be of velocity B 


0.5. Expanding A^(p) to first order in p about we have 
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where represents the velocity derivative of A. (6) evaluated at p = 

We proceed with a two parameter analysis as described in the 
previous section by writing equation 11.16 again for a second observable 
j =f^ i. Solving for the mass, we then obtain the result 


/! = 


m ^ 'Mis: 






(11.17) 


where we have put = p^^ = Ihe uncertainty in the result of the 

two parameter analysis as given by equation 11.17 enters through the 


dependences of the functions A^(p) and Aj(p) on the observations ^ and 
Oj, respectively. In particular, if we let p^ be the velocity correspond- 
ing to A so that A = A^(p^) = A^Cp^), then the uncertainty, in the 

mass determination, A, is 'given by 



(11.18) 


(11.19) 


defines the mass resolution scale factor, and cr^ denotes the resolution 
of the observation In order to interpret the scale factors, assume 

both parameters contribute equally to the resolution in equation 11.18. 
In this case, the condition for the existence of a valley between two 
neighboring mass distributions is given by 


£ 






( 11 . 20 ) 
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s 

where is defined by equation 11.12. In the particular case where 

the abundances of the two isotopes are in the ratio 2:1 (i.e. k = 2) , 
this reduces to 


(n < o,z^ 




( 11 . 21 ) 


max 

The are listed in table II. 2 for the various possible observables 

10 

and analyses, again for the particular case of a Be with velocity g = 
0,5. We also recall that this table strictly applies only to the case 
where k = ^Be/^^Be = 2, and both observables in a two parameter analysis 
contribute equally to the uncertainty in the resultant mass, 

C . The Experiment 

1, Specialization to ExR . We chose to do an experiment measur- 
dE 

ing the observables E, “ and R which utilized the existing technology 

OX 

in large area scintillation counters at our laboratory. The measurement 
of range was to be the unique feature of the experiment. This was accom- 
plished by stopping particles in a stack of relatively thin scintillation 
counters, which then provided multiple measurements of both and E as 
functions of the residual range of the particle. The result is a redun- 
dancy of mass determinations, which may then be used for the identifi- 
cation and removal of background. In the remainder of this section, we 

will specialize the discussion of section B to an experiment of the 
dE 

multiple X E X R variety. 

In figure II. 9, the energy of a proton is given directly as a 

function of its range. Curves corresponding to other particle species 

will have the same (log-log) shape, onl> offset in both R and E according 

to the scaling laws given by equations II. 1 and II. 4. It is apparent 

-2 

from the figure that, over a considerable range (up to R S£ 100 g cm ), 
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TABLE II. 2 


Resolution Required to Resolve 

a Be 




of 



Velocity p 

=0,5 in a 

Two Parameter Analysis 




Primary Observable 


Auxiliary 


dE 



Observable 

E 

dx 

R 

P 

E 

— 

1.77o 

1.77. 

1.37o 

iE. 

dx 

2.37, 

— 

2.37c ■ 

2.57. 

R 

1.07. 

1.07. 

— 

1.87. 

P 

2.37, 

3.37. 

5.87. 




E (MeV) 
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the relation 




= 3Jl,3 



( 11 . 22 ) 


is an adequate approximation to the range energy relation. Solving for 


the mass. A, in terms of R, E and Z, one has the result 




r r.U‘' 


(11.23) 


The resolution of the mass determination is then given by 



(11.24) 



(11.25) 


for the particular case of a Be with velocity p = 0.5. This is in 

agreement with the result given by equation 11.18. Equation 11.24, 

however, applies to particles of arbitrary Z, A and ft. It is apparent- 

from this result that an experiment of the E x R variety has the poorest 

resolution for the heaviest isotope at the lowest energy. 

2. Range Resolution . From equation 11.24, we estimate that the 

range must be resolved to within ~ 27 c in order to clearly resolve the 

isotope (as defined by equation 11.21). Thus if we take the 

_2 

resolution to be o’,, ~ 1 g cm (i.e. the thickness of a thin scintillator), 

K 

1 n ^2 

we conclude that Be can only be resolved for R ^ 50 g cm . Clearly, 

we must do better by at least an order of magnitude if this approach is 

to be feasible (e.g. the mean free path of ^^Be against nuclear inter- 
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actions in scintillator is ~ 22 g cm ), Thus, we are required to 
resolve the range to tnuch less than the thickness of a single scintil- 
lator. A simple solution is to use the- known range energy relation 
(i.e, see figure II. 9) to deduce the range in the last detector pene- 
trated from the energy deposited. If we approximate the range energy 
relation by equation 11.22, then 


['Aj 




l.7f 


(11.26) 


where is the- pulse height (in MeV) in the last detector penetrated. 

This, however, requires an a priori knowledge of A and Z. This problem 
may be circumvented by writing equation 11,26 again for the range in 
the last two detectors penetrated: 

Here thickness of the next to last detector penetrated, 

and 9 is the zenith angle at which the particle is incident. Equations 
11.26 and 11.27 may then be solved to give the result 



ILt-/ 




(11.28) 


which depends only on the pulse heights in the last two detectors, the 
exponent of the range energy relation, the thickness of the next to last 
detector and the zenith angle, all of which are observable. The problem 
is, then, apparently solved, as the range resolution has been consider- 
ably improved. 
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3. Scintillator Saturation . There is yet another effect which 
must be considered, namely, the saturation of the response of the scintil- 
lation counters used. This has two effects. In the first place, the 
response of a scintillator is now a nonlinear function of the energy 
deposited, as well as Z and A. Although the effect is not completely 

understood in a quantitative sense, it does seem possible to conclude 

dL 

that the saturation of scintillator response (i.e. depends only on 

Z and p, and not on mass. The detailed form of this dependence, however, 

must be measured for each specific application. This may be done by 

using an accelerator for extensive calibrations. Alternatively, it may 

be accomplished (with less precision) by using in flight data. The 

latter possibility relies on particles known to stop near detector 

boundaries, for which the range is well known. The second effect of the 

scintillator saturation concerns our ability to estimate the range in the 

last detector. One may hope that an assumption of a power law response 

is still approximately valid, in which case equation 11.22 (with adjusted 

parameters) may still be applied, perhaps in an iterative fashion. There 

is some hope for this since the saturation is expected to be a relatively 

small effect for the particles of interest (i.e. particles of low charge, 

10 

and in particular. Be). However, the exact effect of scintillator 
saturation is not known a priori, and the calibration may be required. 

4. Additional Considerations . Yet another factor which should 
be considered is the requirement to obtain a statistically' significant 
result. As an example, consider an isotope such as ^^Be. The differ- 
ential flux of this isotope is expected to be ~ 5 x 10 ^ particles 
m“^ster“^sec“^ (MeV/nuc) ^ (Tsao et al. 1973). Assuming an exposure 
time of 24 hours and an energy window of ~ 400 MeV/nuc in which the 
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experiment is sensitive, we expect to observe roughly 1700 G particles, 
where G is the geometrical factor of the experiment in ster. Even 

O 

with a geometrical factor G - 0,25 m ster (which is the largest we can 
attain using the available gondola), we expect only ~ 425 ^*^Be, even 
before considering losses due to nuclear interactions, dead time and so 
on. Clearly these effects will be important, so that we should make 
every effort to maximize the geometrical factor. Finally, a constraint 
on the upper limit was set by the requirement to utilize existing hard- 
ware wherever possible, including certain detectors, pulse height 
analyzers and the gondola, 

5, Description of Experiment . To summarize the preceding dis- 

dE 

cussion, the experiment was designed to operate in a multiple ^ x E x R 
mode in order to maximize its isotopic separation. The heart of the 
experiment consists of a stack of scintillation counters with which it is 
possible to determine the range of a particle and also to test for back- 
ground by consistency checks. 

a . Overview 

The physical layout of the experiment is illustrated 
schematically in figure II. 10. Briefly, it consists of a thin trigger 
scintillator, Dl, followed immediately by a Cerenkov radiator, D2. This 
is followed by a spark chamber, which was used to define particle tra- 
jectories through the experiment. This was especially important for the 
identification of stopping particles, the determination of particle 
ranges from their normal projections, and mapping of positioixal variations 
in the responses of the various detectors. The spark chamber is followed 
by the totally active stack of 12 scintillation counters, D3-14, which 
comprised the heart of the experiment. These allowed multiple measure- 
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Figure II. 10. Schematic diagram of the experiment. 
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ments of E and as functions of R for each event analyzed. The final 
element, Dl5, is a penetration counter, and gives only a "yes" or "no" 
output. The scintillators are all made of Pilot Y, chosen for its good 
response to heavily ionizing particles. 

b. Testing and Checkout of Detectors 

Each detector assembly (consisting of scintillator, light 
pipes, photomultiplier tubes and pre-amplifiers) was evaluated in a 

230 

light tight box using muons and a l/iC, P alpha source. For the 

h 

scintillators, commercial amplifiers, pulse height analyzers and high 
voltage supplies were used for these tests; in general, the high voltage 
was set to levels well above the flight values. For the Cerenkov 
detector, the flight amplifier and high voltage supplies were also 
incorporated in the checkout; after balancing, the photomultiplier high 
voltage levels were fixed and not changed for the flight. The general 
procedure of the evaluation was first to infer the response per photo- 
electron of the system from the location, V , and full width at half 

a 

maximum, W^, of the peak from the alpha source, assuming purely Poisson 
statistics. The most probable number of photoelectrons produced for 
each muon, Npgj was then deduced from the location of the muon peak, 

V^, relative to the alpha peak: 


I 



(II .'29) 


The procedure was similar for the Cerenkov detector, except that a 

n / -1 

O.OOSfiC, alpha source embedded in a Nal scintillator was used as 

a light source. The general results of the preflight checkout are given 
in table II. 3, together with other miscellaneous data pertaining to the 


various detectors. 
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c. Detailed Description of Detectors 

The Dl scintillator was included primarily for triggering 
purposes. Its size was made significantly larger than the other scintil- 
lators in order to increase the effective geometry of the scintillator- 
Cerenkov telescope (Dl-3) for measuring the charge composition of non- 
stopping particles. Because of this large size, it was occasionally 
possible to trigger the experiment with particles which did not pass 
through the spark chamber, which was only 50,8 cm x 50.8 cm. Such events 
were discarded in the analysis.- The Dl scintillator was made as thin as 
possible (0.635 cm) in order to minimize the amount of matter above the 
totally active stack of scintillators, D3-14. This was balanced with the 
necessity to have a reasonable level of resolution (dominated by Landau 
fluctuations) for triggering the experiment. 

The Cerenkov radiator, D2, is a sheet of acrylic plastic 
1/2" thick. It is enclosed in a box painted with diffuse reflecting 
white paint, and viewed through the sides of the box. 

The Cerenkov detector is followed by a digitized ■ spark chamber 
consisting of eight decks 1.27 cm apart (Ehrmann et al., 1967). Each deck 
consisted of one x and one y plane of 200 wires each. The wires, spaced 
0.254 cm apart, where threaded through magnetic cores whose' polarity 
could be reversed by a current pulse in the wire. The set cores were 
read out using the pulse induced in a second wire when the cores were 
reset through a third wire. The spark chamber was used to determine the 
pathlength of a particle in a given detector, to correct the data for 
positional variations in the response and thickness of the detectors, 
and also to predict the point at which a particle would exit from the . 
•totally active stack, D3-14, assuming it did not stop. The latter was 
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essential in detertnining whether or not a given particle had actually 
stopped in the experiment. The spark chamber was also useful in reject- 
ing background, which appeared as multiple tracks (Arens et al. 1974); 
events which gave no track were also discarded. Both the spark chamber 
windows and the Cerenkov light diffusion box were constructed of aluminum 
sheets 0.124 g/cm^ thick. 

The D3 scintillator is primarily a trigger scintillator. How- 
ever, it also serves as the first element of the totally active stack of 
scintillators D3-14, and is especially useful for particles of low energy 
which do not penetrate very far into the stack. Its size is the same as 
the active area of the spark chamber, and thus is somewhat smaller than 
the other scintillators. This helped the experiment to trigger on 
particles which had actually passed through the spark chamber. D3 was 
made fairly thin in order to maximize the number of detectors penetrated 
by low energy particles before stopping, thus maximizing the information 
available for the analysis. The limiting factor is again the role of 
statistical fluctuations, together with the physical requirement that 
the experiment should not have large gaps of empty space between detectors 
(which would drastically decrease the geometric factor). 

The three thin scintillators, D4-6, are slightly thicker than 
D3 in order to minimize empty space between detectors; However, they are 
still sufficiently thin to give multiple data points for low energy 
particles. D4r6 have dimensions somewhat larger than D3 in order to 
increase the geometric factor as much as possible within the limitations 
imposed by the size of the spark chamber, which is fixed. 

The thick scintillators, D7-14, are intended to extend the 
sensitivity of the instrument^ to higher energy particles. Although 
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detectors D7-14 are much thicker than the other scintillators, they are 
.still sufficiently thin to provide significant information in the critical 
region near the end of a particle's range where ^ is rapidly changing. 

The anticoincidence scintillator, D15, is .included to give 
some indication of whether or not a particle has stopped in D14; without 
it, D14 is reduced to an anticoincidence scintillator for D13, Particles 
were defined to have stopped in D14 if, and only if, their trajectories 
passed through D15 and D15 gave no response. 

,d. Light Pipes 

There are two types of light pipe used in the experiment. 
Adiabatic,, plastic finger light pipes were used with Dl and D3 in order 
to maximize light collection efficiency for these two very thin scintil- 
lators. For D4-14, triangular light pipes were used. Their shapes were 
determined by connecting the edge face of a given scintillator to a 
rectangle 7.62 cm high (determined by the diameters of the photomultiplier 
faces) by four planar surfaces. They were made as long as, possible 
(within the constraint that they fit into the gondola) in order to max- 
imize their efficiency. They were of two general types, as described .in 
table II. 3. Each type was compared with the adiabatic light pipes and 
found to have about half the collection efficiency of the adiabatic 
pipes. On , the other hand they were much less costly to construct. The 
main problem connected with the use of these light pipes was a rapid 
variation in collection efficiency near the junction of the two pipes, 
which was .corrected in the data reduction. These light pipes also 
resulted in some problems during the balancing of the photomultiplier 
high voltages due to geometrical as 3 nnmetries present in the design. 

e.. Coincidence Requirements 
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A schematic of the logic used in triggering the experiment is 
shown in figure II. 11. As can be seen, the experiment was designed to 
trigger in four separate modes: the Isotope mode (for stopping particles), 

the Cerenkov mode (for heavy particles penetrating Dl-3) , the penetrating 
mode (for events to be used in mapping positional and temporal variations 
in detector response) and the calibration mode (for use in normalizing 
the responses of the various detectors). The criteria required for 
triggering in each of these modes are given in table II. 4. 

Most events triggered the experiment in the isotope mode. This 
was the primary triggering mode for the experiment, and was intended to 
detect particles heavier than helium which stopped in the experiment. 

Since ionization energy loss decreases with increasing range, this mode 
was designed to require lower thresholds on Dl for particles of longer 
ranges. These thresholds were set to reject He (which would otherwise 
have swamped the experiment) while accepting as much Li as possible, 
with no rejection of Be. 

Of the other inodes, the Cerenkov mode was intended as a separate - 
experiment to measure the charge composition using detectors Dl-3, while 
the penetrating and calibration modes were intended to provide calibration 
data for calibrating the various data and removing systematic errors. The 
calibration mode failed to operate properly in the experiment, and so 
provided no useful information, 

f. Pulse Height Analysis 

There were two types of pulse height analyzers used in the 
experiment. All used automatic gain switching in order to obtain the 
required dynamic range of ~ 10^. Detectors D4-14 each had two 256 chan- 
nel analyzers with synchronous clocks and three gain ranges each. The 
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Figure II. 11. Schematic illustration of triggering 
logic used for the experiment. 
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gain, ranges of the two pulse height analyzers *^n a given detector were 

overlapping, thus permitting an accurate calibration of gain factors 

and zero offset from in flight data. Dl-3 used 1024 channel analyzers 

with as 3 nachronous clocks and four gain ranges each. D3 also was analyzed 

by a 256 channel analyzer, which was included primarily for' the extra 

thresholds it provided for triggering purposes; however, it also served 

as a partial check on the 1024 channel analyzer. A schematic diagram 

of the pulse height analysis strategy is given in figure 11.12, 

Flight . The experiment was flown by balloon from Thompson, 

Canada, on August 15, 1973, at 5:04 a .m,. local time. Data was received 

at Thompson and Ft. McMurray. The payload floated at altitudes between 
-2 

3.5 and 5 g cm residual atmosphere for ~ 22 hours and was cut down at 
7:00 a.m. on August 16. 



2 BITS 
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Figure 11.12. Pulse height analysis strategy for detectors D7-14, 







CHAPTER III 


- DATA- ANALYSIS 
A. Formulation of Response 

The differential response of detector i to a particle of charge 
Z, mass A and velocity p is given by the relation 
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where the factor M describes the variation in response with position 
-»• , dL 

<r^) and time (t) for detector i, and ^ represents- the efficiency with 
which energy lost by the particle in detector i is converted to signal. 
Here it is assumed that the systematic variations in response with time 
and position in a detector (represented by may be separated from the 
dependence on parameters (i.e. Z, A and p) which characterize the particle 
itself. We will also assume that ~ is independent of mass and detector 
number. Although this latter assumption (namely, independence of 
detector number) is not justified on a broad scale, it is a satisfactory 
approximation for the particles of charge with which we will be 

exclusively concerned in this thesis, as will be shown later. Finally, 
since the detectors are in general quite thin with respect to their 
horizontal dimensions, any dependence of on r-^ • (and consequently 
also on the zenith and azimuthal angles, 0 and 0, respectively) can be 
neglected. If is the particle’s range (as measured from the top face 
of detector i) , then the pulse height measured by the detector is 
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wtiere T.(r.) s: T.(.lr-) gives the detector thickness as a function of 
i i i j- 
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position, j^^i» in its central horizontal plane. For a relativistic 
particle, this becomes 
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MxTi (III. 3) 

which is the basic relation used for calibration and mapping of position- 
al and temporal variations. 

The response of the Cerenkov detector may be characterized in a 
similar fashion. The differential response is given by 



(III.4) 


with 






(III. 5) 


for a Cerenkov radiator with index of refraction, n . Analogous to the 
case for scintillators, we have 
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(III. 6) 


and 




(III. 7) 


The pulse heights, are of course statistical quantities. As such 
they are characterized by distribution functions rather than unique 
values. These distribution functions are determined by both physical 
and instrumental factors. For example, important contributions come 
from Landau fluctuations in the production of delta rays and Poisson 
fluctuations in the conversion of scintillator light to photoelectrons 
in the photomultiplier. All distribution functions are (for convenience, 



and unless otherwise noted) approximated by Gaussian distribution with 
widths characterized by their standard deviation, a. Furthermore, these 
distributions are all assumed to be statistically independent (again, 
unless otherwise- noted) . ■ 

The pulse height represents the signal (light) observed in 

detector i by a particle of charge Z, mass A and velocity p at the top 

of the detector. Because of the large dynamic range required (typically 
5 

~ 10 ) , this was analyzed by a pulse height analyzer with dynamic gain 
switching. There were two types of pulse height analyzer used: D4-14 

were each analyzed by two 256 channel- analyzers using three gain ranges 
each, while Dl-3 were analyzed by 1024 channel analyzers using four gain 
ranges each; D3 was also analyzed by a 256 channel analyzer. Since each 
gain range can be regarded as a separate pulse height analyzer with its 
own gain and zero offset,- the pulse heights are related to the raw 
channel numbers, H'-, (where j denotes the gain range) according to 



(III. 8) 


The determination of the and bj^j is discussed in detail in Appendix A 
and so will only be briefly sketched here. The and were all 
determined by terms of and b^^^ by use of the overlapping gain ranges 

of the two different pulse height analyzers analyzing a given detector. 

In cases where there was only one pulse height analyzer, b^j = 0 was 
assumed, and the were determined in terms of by comparing the 
cut-out channel of a given gain range with the cut-in channel of the next 
gain range (these should both correspond'to the same pulse height). 
Alternatively, when two pulse height analyzers analyzed the same detector 
both the a^- and b^- could be determined in terms of the a^-j^ and b^j^ by 
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plotting one pulse height analyzer against the other. Next, the and 
bii were all determined in terms of and using semirelativistic 
particles (i.e, particles incapable of producing knock-on electrons of 
long range which might couple one detector to another, yet sufficiently 
fast that they did not slow down significantly in traversing a detector). 
Finally, b^j^ = 0 was assumed, and a^j^ was determined to give reasonable 
pulse heights for relativistic particles (i.e. in x •minimum units). 

In addition to the pulse height analyzer gain and zero offset 
factors (a., and b^O, the factor H. (which describes positional and 
temporal variations in response) must also be determined. M^(rj|^,t) has 
been assumed to be separable: 

,*) = flfd) 

P 

Here the mapping function, M^, depends only on position and the drift 

T P T 

function, depends only on time. Both and are normalized to 

P 

unity in order to preserve detector normalizations. The functions 
T 

and were determined by observing variations in detector response to 

P 

relativistic C and 0 nuclei. was approximated by an 11x11 polynomial 

T 

with the 121 coefficients determined by a maximum likelihood fit. M., 
on the other hand, was approximated by the least squares spline technique 
(Thompson 1973) . 

Finally it was necessary to know the thicknesses of the various 
detectors as a function of position. This was accomplished by caliper 
measurements of each detector on a 5 cm grid. 

B. Preliminary Selections 

In the preceding section, the reduction of the raw data to a form 
essentially independent of the idiosyncracies of this particular experiment 
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has been discussed. The data must new be analyzed in a manner capable of 
achieving elemental and Isotopic resolution. This has been approached 
in different ways for various subsets of the data: For very high energy 

events (E Js 570 MeV/nuc), only the elemental abundances have been derived 
using a Cerenkov versus scintillator analysis based on detectors Dl-3. 

The elemental abundances resulting from this analysis are shown in 
figure III.l, Events which stop in D3-14, and yet are above the Cerenkov 
counter threshold (i.e. 5 ^ 0,67) are subjected to an isotopic analysis 
based on the Cerenkov and range measurements. This has the unique 
advantage of not requiring the determination of any peculiar response 
functions — the interdependence of Cerenkov response and the range of a 
particle are well known. This approach is particularly effective for the 
heavier particles with Z ^ 8 . Finally, particles which stop in D3-14 and 
have 4^Z^7 have been subjected to a light versus range analysis. The 
restriction on the charge of .the particles analyzed is imposed by trig- 
gering thresholds from below, and from above by the non-uniform response 
functions of the Pilot Y scintillators used in the experiment. This 
thesis will be concerned exclusively with this mode of analysis, which is 
described in detail in section B of the present chapter. 

The data analysis begins with the isolation- of a relevant subset 
of the data, which is defined by the following criteria. 

Only events stopping in detectors D5-13 are included in the 
analysis. A particle is said to have stopped in a given detector if that 
detector produces the last nonzero pulse height observed and the extra- 
polated spark chamber track exits from the stack at a point deeper than 
halfway through the following detector. Events stopping before D5 suffer 
from reduced resolution due to the minimal number 'of detectors penetrated. 
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Dl4, on the other hand, suffers from noise from the spark chamber in one 
of its pulse height analyzers. In addition, it is not always possible to 
determine whether a particle has stopped in this detector due to the 
reduced size of the anticoincidence scintillator, D15. 


No events producing either no track or more than one track in the 
spark chamber are included. An acceptable track must include discharges 
in at least five of the eight spark chamber decks in the x and y views. 

"Obviously interacting" events were excluded. These were 
identified by examining the sequence of pulse heights in detectors D3-14. 
Events with pulse heights more than 157» less than the average of the 
pulse heights from the preceding two detectors (excluding the Cerenkov, 
D2) or more than 30% less than the pulse height from the immediately 
preceding detector were excluded unless the pulse height in question 
corresponded to the detector in which the particle stopped. 

Finally, events were discarded in which the pulse heights from 
two pulse height analyzers analyzing the same detector did not agree to 
within one channel equivalent of the smaller of the two pulse heights. 

C. R X L Mode of Analysis 

The R X L mode of analysis (with which we are exclusively con- 
cerned here) is based upon a presentation of the response of the totally 


active stack of scintillators in terms of the variables L£(total light 
observed in and after detector i) and (range of particle measured 


from top of detector i) . In doing this, the implicit assumption that 

the detectors all have identical response functions, has been made 

d£ 

in order to allow the interpretation of as a simple sum of pulse 
heights in different detectors. As has been observed previously, this 


assumption is noc generally valid, and is a suitable approximation only 
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over a limited range of incident particles such as that (4^^7) with 
which the present analysis is exclusively concerned. As an example, the 
responses of the various detectors to relativistic particles of various 
charges are tabulated in table III.l. Ideally each detector should show 
the same response to the same incident particle in this table. As can 
easily be seen, however, this is not the case, and indeed cannot be the 
case over the entire range of charges for any choice of detector normal- 
izations. However, for the restricted range of particles considered 
here, the effect is 17i. The R x L mode of analysis offers important 
advantages due to the scaling of the fundamental parameters and R^ 
in mass. It is an intuitively straight forward approach, which still 
offers the advantages of a multiparameter analysis for testing the 
Internal consistancy of events in order to identify and eliminate back- 
ground. The R X L analysis is very similar to the R x E analysis which 

was discussed in Chapter II. Its fundamental concepts are illustrated 

9 

in figure III. 2, which depicts its application to an actual Be event. 

In figure III.2A, the event is shown incident on the experiment at a 
zenith angle of 29°, and stopping in detector Dll. The analysis of this 
event will now be described in detail in the following subsections. 

1, Estimation of Range . As is shown in the lower half of 
figure III.2A, one first adds the various pulse heights to obtain several 
values of (one for each detector penetrated) as a function of detector 
number. Insofar as detector numbers may be related to depth in the stack, 
this is already a crude R x L representation. As we have seen in 
Chapter II, however, this estimate must be considerably refined if we 
are to have isotopic resolution. The refinement of the range estimate 
obviously depends upon a determination of the range in the last detector 



TABLE III.l 


Responses of Detectors DI and D3-*ll to Relativistic 
Particles of Various Charges ■ (xMlnimum Units) 
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Figure III. 2. Schctnacic illustrntlon of range versus 
total light analysis. 
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(in which the particle stops), Even' without such an estimate, 

however, there is a certain subclass of particles for which the range is 

fairly well known, namely those which stop near detector boundaries. 

Such particles may be identified by looking for (relatively) very small 

pulse heights in the last detector’ penetrated by stopping events, i.e. 

^last"^ 0. For such events (depending upon the limits put on 

R may be made arbitrarily small and, consequently,, negligible compared 
Xd s u 

to detector thicknesses (which are added to it in the analysis). Of 
course, the tightness of this selection must be balanced with the 
statistical weight of the resultant subset of the data. Figure III. 3 
shows the result of plotting such a subset of the data on what amounts 
to a rotated In L versus In R plot (i.e. the dependent and independent 
variables have been chosen as linear combinations of In L and In R such 
that the response curves are nearly horizontal). As can be seen, the 
detector response is well approximated by a power law of the form 




(III. 10) 


where 6 1.3 is a constant, and qt(Z,A) does not depend on velocity. By 
writing this equation again for the next to last detector penetrated, the 
parameter c^(Z,A) may be eliminated. One may then use the relation 


W-1 






(III. 11) 


to obtain the result 




Vst 


= T 3.6 


(III. 12) 


from which it is possible to estimate the range of a particle in the last 
detector with no a priori knowledge of its charge or mass. The derivation 
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Figure III. 3, Rotated light versus range plot for particles stopping 
near detector boundaries. 
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of this relation requires only that the (integral) response of the last 
two detectors to a given particle be characterized by identical power 
laws of the form III. 10. This assumption appears to be justified at 
least for 4^:e7 and 






> 

rJ 




0 ,^ 1 ? 3 



(III. 13) 


where easy experimental checks are possible (see figure III. 3). For 
particles of shorter range, this dependence must then be extrapolated. 

From a strictly mathematical point of view, it is possible to 
extract some information on the detector response, even for R < T^. 
Essentially, this would be done by requiring 








(III. 14) 




T seed ((- 




— *r sgc. 0 


which is based entirely on observables. This would be useful as a cross 
check on the reliability of the range estimation procedure. However, both 
terms in III. 14 are generally larger than their difference, with the result 
that the accuracy suffers. For this reason, and statistical limitations 
imposed by the amount of data available, it is not possible to use 
equation III. 14 to meaningfully extend the checks on equation III. 12. 

An attempt to compare the predictions of equations III. 14 and III. 12 is 
shown in figure III, 4. 

2. Standard Curves . Having determined first order, all 

of the available data may be presented on a diagram similar to that of 
figure III. 3. This is shown in figure III. 5, From this figure, one sees 
that there is a very clear separation of the various charges. However, 













/nL -o.fr/nH 



Jol + /nR/0.77 


Figure III. 5. Rotated light versus range plot for unselected data with standard 
response curves sketched in. 
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mass separation is apparent only for Z ^ 3 and Be. In particular, the 

10 9 

important isotope, Be, is not clearly separated from Be, and the 

problem gets worse for heavier isotopes. Since the response function 
dL 

is independent of the particle mass (A) while depending on both charge 
and velocity, it should be sufficient to determine ^ for a single isotope 
of each element, and then extend the result to the remaining isotopes. 

This will be discussed later in more detail. 

The standard response functions L(R) are obtained by averaging 
over the widths of the .distributions in figure III<,5 along the ordinate. 
Where at least one isotope of a charge group is clearly separated, only 
data pertaining to the most abundant such isotope is used. Where no 
isotope is clearly separated (i.e. for Z > 5) all of the data for a given 
element are treated as if pertaining to a single isotope. This does 
carry with it certain difficulties. For instance, the density of experi- 
mental points along a given response line in figure III. 5 is a compli- 
cated function of spectral characteristics and the geometry of the 
telescope. In reality, the vertical contours used to locate the standard 
response curves are probably not always normal to the gradient in the 
density of experimental points. It is thus probable that there are 
small systematic errors in the determination of the standard response 
curves for elements with significant abundances of more than one unresolved 
isotope. The effect is most important for B and N (which do not exhibit 
a single, dominant isotope), and does not apply to Be at all since the 
Be response curve is based on the well separated 'Be data. Presumably, 
it can be corrected by an iterative approach if the initial approxi- 
mations to the response lines are sufficiently close to reality. However, 
this has not been attempted in the present analysis. 
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The standard response curves which have been adopted are shown 
in figure III. 5, together with the nominal masses assigned to them. The 
precise masses to which they correspond will be determined later in the 
analysis. In the particular case of B, the response curve has been 
shifted somewhat (according to rules to be developed in the following 
section) to correspond more closely to the isotope These same 

standard response curves are presented in figure III. 6 in the form of 
integral response (saturation) as a function of velocity by use of -the 
range energy tables. Also shown are the corresponding differential 
response functions. 

3. Mass De terminat ion . In order to apply the standard response 
lines to the data, use is made of two fundamental scaling laws. The 
first. 




(III. 15) 


is a restatement of equation II. 4, while the second, 




(III. 16) 


is analogous to equation II. 1, Its validity is equivalent to the assump- 

dL 

tion that the differential scintillator response function, , is independ- 
ent of the mass of the ionizing particle (see equation II. 3). This is 
reasonable, at least to first order, since the medium sees only the 
charge of the incident particle. For Coulomb collisions with atomic 
electrons, the mass of the projectile (assumed to be an atomic nucleus) 
is relatively very large and does not significantly affect the energy 
transferred to the knock-on electrons or the number of collisions per 
unit pathLength which are the two fundamental parameters involved in the 



Figure III. 6. Integral and differential saturation as 
functions of velocity as derived from 
standard response curves. 
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problem (Rossi 1952). The most plausible objection to this scaling would 
be to note that saturation (i.e. the departure, from unity of the scintil- 

ik 

lator response, is a macroscopic phenomenon (which is not well under- 

stood) rather than a microscopic one (Meyer et al. 1962). The suggestion 
is then that one should take account of the fact that a particle of smaller 
mass slows down more in traversing a macroscopic distance (~ the dimension 
of a "typical luminescense center") than one of greater mass and identical 
initial velocity. Thus the lighter particle will deposit slightly more 
energy and the scintillator response may be more highly saturated. 
Essentially, the condition for this effect to be negligible and the 
scaling given by equation III. 16 to hold is 




« 




d itni -bI 

d. iL/n R 




(III. 17) 


where the second inequality results from the power law approximation to 
the range energy relation (equation 11.22) and the response curves given 
in figure III. 6. As an upper limit, we may take )Li to be a character- 
istic dimension for luminescence centers. A better upper limit, fx' , would 
be the magnitude of the range discrepancy which develops when two particles, 
of the same charge and velocity but differing mass, traverse a distance 
equivalent to the characteristic' dimension of a luminescence center, /i: 



(III. 18) 


Here and ^2 given by 



/^Z 

I 




(III. 19) 
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and is the initial velocity o£ both particles. In any case, the 
inequality III. 17 is easily satisfied for any realistic values of the 
relevant parameters. A more detailed experimental and theoretical invest- 
igation of the phenomenon of scintillator saturation has been given by 
Taylor et al. (1951), Meyer et al. (1962), Katz et al. (1968), Kobetich 
et al., (1968). 

The analysis now proceeds to figure III.2B. Using' the scaling 
law given by equations III. 15 and III. 16, we obtain 

Similarly, 

(III. 21) 

Thus, on the log log plot of figure III.2B, the curve corresponding to 

variable mass with fixed charge and velocity is a straight line of slope 

1. Furthermore, the distance between an experimental point (mass A) and 

the standard response line (characterized by mass Ag) along such a (slope 

A 

1) curve is just yiln a . In this way, a mass can be defined for each 

experimental point. This, of course, depends upon a priori knowledge of 

the charge of the particle; however, as is apparent from figure III. 5, 

the charge determination is straight forward. Should the charge not be 

well defined, one would simply compute a mass for each possible charge, 

then select the charge for which the various mass estimates are most 

consistent. Finally, since the points ^last^ ^^last-1’ 

R, -) have already been used to define the residual range in the last 
last-1 

detector, they can provide only one unique estimate of the mass. The 
results of the analysis just described for the event considered in figure 

III. 2 are shown in figure III.2C, 
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The best estimate of the mass is taken as that corresponding to 

the maximum range, namely This is reasonable since, in general, the 

point (L^,R^) contains all the information contained in pl^s 

a contribution from detector i. As one goes to longer and longer ranges, 

the increase in precision due to additional detectors penetrated decreases. 

In other words, the most important detectors in determining the mass of a 

particle are the ones nearest the end o-f its range. This especially 

dL 

stands out when one chooses to analyze in a multiple ^ versus (R or L) 

dXr 

modei where 35^ must be measured as close as possible to the end of a 
particle's range. The mass estimates A^, , . are used only to 

test the internal consistency of the event or, in other words, to look 
for background. 

4. Rejection of Background . One would expect the to approach 
A3 in the manner of a convergent series as detector number, i, decreases. 
This follows since, as range is increased, each new detector added con- 
tributes a smaller fraction of the total light, L^. Such a behavior is 
seen, for example, in -figure 111.2c. If we suppose that the are dis- 
tributed about A^ with a Gaussian distribution, then we can compute the 

2 

variance of the A^_, var(A) = cr , where ct is a number characteristic of 
the width of the distribution. However, for background events, one 
expects relatively large excursions of the since the particle has 
changed its mass (and possibly also its charge) at some point in the 
experiment. These excursions should be reflected in abnormally large 
values of the var(A) parameter. In figure III. 7, a frequency distri- 
bution of the var(A) parameter for Be events is shown. From this figure 
a test given by 
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Figure III, 7. Distribution of mass variance for Be events. 
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has been defined in order to select non-interacting events. Events which 
fail this test are identified as background and discarded. The scaling 
with mass of the selection threshold given by equation III. 22 can be sup- 
ported on both analytical and observational groundsj this justification 
will be discussed in detail when the mass resolution is considered in 
section III.D. The effectiveness of the criterion ill. 22 in removing 
background events is readily seen by a comparison of figures III. 8 and 
IIIo5. 

Observed Isotopic Composition . The mass histograms resulting 
from the analysis just described are displayed in figure III. 9. Also 
shown are equivalent histograms for events which were rejected by the 
var(A^ criterion. III. 22, as background. It is to be noted that, in most 
cases, the background level is quite flat; for the one notable exception 
(C), the background peak is at a distinctly different mass than the non- 
background peak. As may be seen from the figure, the percentage of events 
rejected decreases as mass increases. One would expect the opposite 
behavior, if anything, since the total inelastic cross section increases 
with A. This observation is taken as an indication of the increased 
efficiency of the initial selections on the data in removing background. 

In particular, the requirement that the scintillator pulse heights be 
(within specified limits and excluding the detector in which a particle 
stops) a monotonically increasing function of detector number is con- 
sidered to be responsible for this effect. 

The mass distribution for each element has been fitted using the 
maximum likelihood method. In doing this, tfie mass distribution cor- 
responding to a single isotope of mass A and charge Z is assumed to be a 
simple Gaussian; 



/nL-0.77 /nR 



5.6 

5.7 
5.6 
5,5 
5.4 
5.3 
5.2 
5.1 
5.0 
4.9 

4.8 


in ! ; .! . •! L i ‘.'ii ,,i i ... ( i 


I X 




j| » MM*I 

i*l» ttn ip »»»*«*o**MM»l****M«*«l«««*a*«**aBai 


IN 




.. . .1 IMUflU a' 




I UM> >MI*M 



... ...^. 
i 1 


> util.fa* ION 


1«M* I ^a.f« > ■■ 

«>a* .all a <a<t.*«a«l«al» aM.I MiaataMtaaaaaa *•<<<■«•«•* | | 

a* 1 1* a*4 ^aaap.lljIM p api •Ma4|«ia|«a.M ppfppappp »«aMV«»iPaM ^ItMaaaai 

M«ira. p H4|*l II aa 44ia»p« raaa.a^ * .ai^a.IaPail apaaaraf n 
4| tt/>*l i ( lallJi <||p.rJtaMa|a»f>aj> *a«aa«a|*4(i|4|4aia4 | | 


& 




I * pa*«a«4 

> ) ( I K»t 


Be: 




• aM'aati ii laiaaaaaii M'aa'iifliatii a I p 

'* * »**ia n III I^A.«a|l .aaaf^aaaaa J| al ii | 

•*■■ - " pPMa'aiaa MfMliiii 


,|aa<a*» .. 

i Ua Maiaiaa Ml IliaapaiaaaPai. ^ aliMt ai 
II <aaaaa>al H PaMii i.Jaaiai Ma^pl^ra a 

IliaJI IIMIkll l^tla I irit»iaiia%Pl(*jl»a» 


t 

P- 

m 






! r 




iai» I 


I M »ti 

i 


As 

tm 


rii 


0^ 

isu 


I 


.. ------- , 4aa«|l^a| 

IMi44ifil«Pta>taMra44a>aai aiaiaai « I i IMM 

a • I aa wra itaiipaMpiaaa will* | | , , , 


a««a» Ma*aiM naa.aaral>«ta«ii 


4 


•r 


a* I aiaaaaa.taaaaaaMPiMiiH 


Ii 


13 


1 


15 


j?n L + jfnR/0,77 


Figure III. 8, Rotated light versus range plot for selected data. 


00 






NUMBER OF EVENTS 


79 




Figure III. 9. Mas$ histograms for data accepted and rejected 
by mass variance criterion. 'Maxiimim likelihood 
fits to accepted data are also indicated. 
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where A is the actual isotopic mass and A' represents the observation. 
Furthermore, the widths of the mass distributions, 2 ;a > different 


isotopes of the same element are assumed to scale as the mass; 


= (A-l ^-*0- (III. 24) 

Z Us) Z'^ 

This is equivalent to the scaling (equation III. 22) of the background 

2 

mass variance threshold with A . The justification for both of these 

assumptions will be discussed later in section III.D. Finally, it will 

be recalled that the analysis only yields the mass of a particle relative 

to that corresponding to the standard response line for a given element, 

A„. It is therefore convenient to include A„ as a free parameter in the 

fit, especially for elements where no isotope was clearly separated. 

The parameters which are varied in the maximum likelihood fit for a 

given element are, then, the mass corresponding to the standard respo.nse 

Ac 

curve, Ag, the characteristic resolution, 2 abundances of the 

various allowed isotopes, best fits to the various distri- 

butions are indicated by smooth curves in figure III. 9. The best fit 
parameters are listed in table III. 2, together with the statistical 
errors (standard deviations) assigned by the maximum likelihood method. 
These error estimates include the effects of correlations among the 


various free parameters in the fit. 

In order to improve the convergence of the maxitmim likelihood 


fit to the mass distribution to nitrogen, a second analysis has been 

14 

performed in which the value of ^cr was fixed. The precise value used 



TABLE III. 2, ■Maxiiaim Likelihood Best Fit Parameters, Uncertainty estimates 

are based on the assumption of Gaussian likelihood distributions. 
The subscript S refers to the standard response curves shown in 
figures III. 4 and III. 5. 
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for ^ has been obtained essentially by generalizing equation III. 24 

to relate resolutions for isotopes even of different elements. Equation 

III. 24 then becomes similar to the assumed scaling of the variance 

selection threshold given by equation III. 22. Both analytical and 

empirical justifications for these assumptions are given in section D, 

where the mass resolution is considered in detail. In particular, 

14 

figure III. 10 has been used to fix the value of The result of this 

procedure, given in table III. 2, is apparently a significant improvement 

in the resolution with which the remaining free parameters are determined. 

However, it must be bom in mind that uncertainties in the fitted para- 

14 

meters resulting from an uncertainty in j through correlations will 
hot be included in the uncertainties estimated by this approach. Never- 
theless, these results have been adopted for the remainder of this thesis. 

In performing the maximum likelihood fits to the mass distributions 
for the elements B, C and N, it was found that they tended to converge 
better when allowed to add contributions from unstable nuclei to the 

mass distributions on the high mass side. The abundances of these isotopes 
12 16 14 

(i.e. B, N and C) were generally more than could ba explained by 
production in the atmosphere above the experiment. The most likely 
cause of such an effect appears to be an asymmetry in the mass distri- 
butions for these elements. This can very probably be related in part 
to the fact that the response curves for these elements were not deter- 
mined from data corresponding to a single isotope. In order to estimate 
the possible magnitude of this effect for Be, we observe that these 
spurious isotopes are generally ~ 10% of neighboring abundant isotopes. 

As applied to this would mean a reduction by one standard deviation 


in the observed abundance. We note that there seems to be no such effect 
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7 8 

in the case of ^Be — in particular, the ~ 19 Be we might expect from the 
above considerations are not observed. For the purposes of further 
analysis of. the data, the spurious isotopes C and have been 

obsorbed into the observed abundances of the remaining isotopes of B, C 
and N, respectively, in such a manner as to preserve the fractional 
isotopic abundances within each element, 

D. Consideration of Mass Resolution 
1, Contributing Factors . The mass resolutions achieved in the 
mass distributions of figure III. 7 are estimated in table III. 2. It 

7 

is readily apparent that only the isotope Be has been clearly separated, 

O 

and this only by virtue of the absence of Be. In particular, we have not 
achieved the goal (o ^ 0.35 for ^^Be) which was set in section III.B.7, 
and consequently do not observe two peaks separated by a valley in the 
mass distribution representing Be and Be. As may be seen from table 
11,2, this may be traced to the uncertainty in our measurement of the 
total light, which is estimated as ~ 2% (we require ~ 1%). This is 
limited by the internormalization of the responses of the various 
detectors, which in turn was limited primarily by the necessity of using 
semi-relativistic particles and also by response functions which varied 
from detector to detector. Thus, the primary goal in improving the 
resolution would be to improve the precision ot the total light measure- 
ment through a detailed calibration of the response of each detector. 
Nevertheless, the resolution which has been obtained is comparable to 
those of other experiments which have been reported (Webber et al. 1973, 
Preszler et al. 1975, Garcia-Munoz et al. 1975, 1975a). 

As a step towards understanding the mass resolution realized in 
the experiment (e.g. see figure III. 9), various contributions are 
illustrated in figure III. 11. The dominant contribution is from Landau 
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Figure III. 11. Contributions to the experimental uncertainties. 
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fluctuations, which result from the statistics of producing very high 
energy knock-on electrons and thereby depositing an abnormally large 
amount of energy in a given detector (e.g, see Rossi 1952). For thin 
detectors, this effect can cause a distribution of pulse heights in 
response to a uniform beam of incident particles which is skewed towards 
high pulse heights. However, for the present considerations, a Gaussian 
approximation to the Landau distribution is adequate. This is given for 
scintillation detectors by 


CP 

LomtloiX. 



1 

J. 







(III. 25) 


where T denotes the detector thickness. Equation III. 25 is valid only in 
the limit where the energy given a single knock-on electron is small in 
comparison to the total energy lost in a detector, which holds for all 
cases of interest here. The magnitude of the Landau fluctuations is 
probably not significantly affected by saturation of the scintillator 
response, since the fluctuations are attributed to long range delta rays 
which deposit most of their energy far from the core of intense ionization 
immediately surrounding the trajectory of the primary particle. However, 
the magnitude of the effect relative to the observed pulse height will, 
in general, increase as scintillator saturation is introduced since the 


pulse height itself is reduced by saturation effects while Landau fluctu- 
ations probably are not. The correction to the Landau fluctuations for 
scintillator saturation effects, derived from figure III. 6 , is also 
indicated in figure III. 11. 

A second contribution to the broadening of the mass resolution 
comes from the Poissonian statistics associated with the fact that only 
a statistical fraction of the photons produced by the scintillators is 
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actually observed by the photomultiplier tubes. This is formulated in 
equation II. 2. The photoelectron conversion efficiency shown in figure 
III. 11 has been taken from the pre flight calibration of the various 
detectors (see table II. 3). 

Of the remaining contributions to the broadening of the mass 
resolution shown in figure III, 11, digitization refers to the round off 
error due to the digitization of the signal by the pulse height analyzers. 
For the 256 channel analyzers (with synchronous clocks), this has been 
approximated by a Gaussian distribution characterized by a standard 
devia tion 


(J“ Ay — ^ (III. 26) 


Estimates of the statistical resolution of the temporal and positional 
mapping of response variations, our knowledge of detector thicknesses and 
the resolution of the spark chamber in defining particle trajectories are 
also Indicated in the figure. Not included in the figure is the contri- 
bution of detector normalization uncertainties, which are thought to 
enter at about the 27. level. 

2. Analytical Analysis of Resolution . After incorporating the 
scaling laws III. 15 and III. 16 into the power law approximation to the 
standard response curves (given by equation III, 10) one obtains a 'general 
formulation of the detector response, which may then be solved for the 
mass : 


_L_ 

/Is 



I 



(III. 27) 


Assuming the independence of the parameters L and R, the mass resolution, 
<J^, may be estimated in terms of the range resolution, and light 



resolution, cr^: 
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(III. 28) 


The assumption that R and L are independent is technically valid only in 
the limit of long range, where the information used to determine R 
becomes an insignificant contribution to L. However, the contribution 
of range to the total uncertainty is normally quite small, as will be 
seen below. It is then reasonable to use equation III. 28 to estimate 
the resolution of the mass measurement even for particles of short range. 

In order to estimate the uncertainty of the range measurement, we 
rewrite equation III. 12 for the range in the last detector in the form 


Rw = Tji.*., -I] 


(III -29) 


where 




Hw-i 


-t I = 


L 




(III. 30) 


gives the ratio of light observed in the last two detectors penetrated 
to that observed in the last detector penetrated. Generally, we have 

2; however when the particle stops in the first thick detector, D7, 
this becomes > 1.25. If the pulse heights are known to a resolution 

O’ II 

given by — then the uncertainty in the range measurement is 
n 



Hr 'i 


(III. 31) 


For H >1.25, the term in brackets is bounded by 1 and 6 and 
R 

equation 111,31 becomes 



(III. 32) 
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Finally, by putting 


a ^ La - K 


JWt 


(III. 33) 


we- obtain. 


^ A. 


S-t I 




(III. 34) 


R SZ[z(Uit-S-) ^ 

The uncertainty in the - determination of the total light, L, is 


given by 




I 


t 


I A ^ 

I H 


(III. 35) 


Substituting the results given by equations III. 34 and III. 35 into 
equation III. 28 and factoring out the contribution from the total light, 
we obtain 


5 

A 






-£-y 


2- 






-■i 


(III. 36) 


The second- term in brackets is negligible for 6^1, which is easily 
satisfied for all cases of interest. In particular, 6 = 1.3 for the 
Pilot Y scintillators used in the present experiment, and 6—1.8 for a 
completely non- saturating scintillator (see equation 11.22)., Thus the 
contribution of range to the total uncertainty in the mass measurement is 
in fact negligible, as has already been pointed out. We then have 


fS-1 

A ^ s '1 / V ^ V H ; 


(111,37) 


For the present case, the scintillators have been estimated to be normal- 


ized to within H 0.02. Putting 6 = 1.3 and assuming the particle to 
H 
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stop in D7, we then have the result 

0.0 31 A CIII.38) 

The result given by equation III. 37 supports the assumption that 
the mass resolution scales as the mass (equation 111,24). It also sup- 
ports the related assumption that the mass variance selection threshold 
used to reject background should scale as mass squared. Essentially, 
the scaling with mass depends only on the assumption of power law response 
functions (equation III. 27), and thus is quite generally applicable. 

The result given by equation III. 37 also gives some idea of the 
parameters one might try to optimize in attempting to design a similar 
experiment with improved mass resolution. For example, the biggest 
improvement (potentially at least a factor two) would result from an 
improved calibration of the detectors. This is readily seen from figure 
III. 11. Such improvements might be accomplished through extensive accel- 
erator calibrations of the various detectors and/or by the construction 
of detectors with more uniform response. Alternatively, one could expect 
as much as a factor of two improvement in mass resolution .by going to 
completely nonsaturating (i.e. & = 1.8) scintillators. To put the 
potential for improvement as just discussed into perspective, we note 
that the mass resolution of the experiment for ^^Be is given by equation 
III. 38 as CT ^ 0.4 amu. This should be compared with the critical reso-- 
lution at which a valley begins to appear between Be and Be. Taking 
k = 2.5 in equation 11.12, this is ‘^crit ~ 0*35 amu. Thus we need achieve 
only a small part of the potential improvement of the mass resolution 
before the results become markedly more convincing. 

Finally, equation III. 36 may be rewritten to give some insight 
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into the problem of deciding the optimal detector thicknesses for this 
type of experiment, namely those such that, the errors in the range 
measurement actually are negligible. In particular, the second term in 
brackets in equation 111,36 represents the ratio of the contributions of 
range and the total light to the total uncertainty. This term may be 
written 



~~Z, 

zizLst -ir/ 


(III. 39) 


where the inequality expresses the condition for the uncertainty in the 
range measurement to be negligible. Putting 6 - l*3j we then require the 
range in the last detector to satisfy 


R., « (III. 40) 

JUit 

Thus, each detector should be much less than half the thickness of the 
material above it in the totally active stack (i.e. D3-14 in the present 
case). This condition is easily satisfied by the present experiment, 
since no particles have been accepted which stop before detector D5. 

3, Empirical Consideration of Resolution . Since the assumption 
of- a mass resolution which scales as the mass is quite important in the 
analysis, it is desirable to obtain empirical evidence of its validity if 
possible. This has been attempted in two ways. The first (see figure 

As 

III. 10) is to plot the ^ a derived from fitting the mass distributions 
against • Since the fit for each element was done independently of the 
others, this gives one point for each element considered. Although the 
resolution appears to scale as predicted, this is not a completely 
definitive check since a scaling of the resolution with charge (rather 
than mass) would satisfy the data equally well. However, there is no a 
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priori reason to believe this is so. A second check was attempted in 

order to test the scaling among different isotopes of a single element. 

Figure III.12A gives the distribution of var(A) for Be with A > 8 as 

compared that for A < 8. If a mass inde pendent threshold on var(A) is 

correct, then the distribution of var(A) should be identical for all A. 

Figure III.12A suggests that this is not the case. Figure III.12B shows 

2 

the same comparison after correcting the variance for the assumed A 
scaling. The situation is definitely improved over the uncorrected case, 
and the mass scaling of the resolution is supported. 

E, Corrections to the Data 

In order to determine the cosmic ray composition at the top of the 
atmosphere, various corrections have been applied. These factors are 
listed in table 111,3, together with the observed and various corrected 
abundances. The corrections treat losses due to interactions in the 
experiment itself (including the gondola), the production and loss of the 
various species in the atmosphere, and energy normalization factors to 
account for the fact that different isotopes are observed over different 
energy windows and modulated differently in penetrating the solar cavity, 
1. Interactions . The correction factors for interactions in the 
detector itself (column C of table 111,3) are given by 




C‘ 



‘iR." ^ (At 


(III. 41) 


Here n(Aj.^j.) is the density of target atoms of mass A^^^ along the 
trajectory of a particle, and the total inelastic (reaction) cross section, 
, is given by the empirical relation 




(111.42) 
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Figure III. 12. Mass variance distributions (both uncorrected and corrected for 
assumed scaling) for Be events with A > 8 as compared to Be 
events with A < 8. 
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where 

y~i 

t= I— 1.18^ a^p AtflLY')j 


(III. 43) 


This relation has been fitted to the data of Lindstrom et al. (1975) for 

^^0 and beams incident on various targets at the Bevatron, toget- 
her with the data of Renburg et al. (1972) for proton beams. Finally, 
the number of particles stopping between R and R +dR, ^ ^tot j 
computed as a function of charge from the data, with mass dependent 
effects (if any) being averaged out. This is required both by the problem 
of statistics, and the fact that most isotopes are not resolved on an 
event by event basis. 

2. Energy Corrections . Since the experiment is sensitive to the 
various isotopes observed in energy windows of differing width and location, 
the results must be accordingly corrected before meaningful comparisons 
can be made with either theory or other experiments. The ideal way of 
presenting the data would be in the form of differential energy spectra, 
in which case the corrections would be minimized. However, this approach 
is of little advantage for isotopes of very low abundance, in which case 
the statistical errors severly limit the resolution obtainable. The 
alternative is to assume the shape of the spectra, which may then be 
used to normalize the results to some fixed energy. Ideally, the spectra ’ 
used should be experimental, but a complete set does not exist which 
covers all isotopes of interest at the right time (i.e. in August 1973). 

We have therefore adopted the local interplanetary spectra for 1973 which 
have been calculated by Garcia Munoz et al. (1975b). They have assumed 
primary source spectra of the form (E + 400)” which are then propagated 
to the vicinity of the sun using a leaky box model for the interstellar 
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propagation (to be discussed later in Chapter IV) with a mean pathlength 
of Xg 5 g cm* and a density of interstellar matter given by n 1 atom 

cm* . They have then used the numerical technique of Fisk (1971) to 
modulate the spectra, and thus are able to predict the spectra observed 
near the earth. Wherever possible, the calculated spectra have been 
matched to observations by adjusting the available free parameters (mainly 
in the diffusion coefficient used for modulation). The resulting spectra 
may thus be regarded as extrapolations of the available experimental data 
to cover all isotopic species. The correction factors for the various 
energy windows in which the experiment is sensitive have been computed 
mapping the energy windows to the top of the atmosphere, then integrating 
the (predicted) spectra over both the extrapolated energy window and a 
standard energy window (233-494 MeV/nuc) corresponding to that in which 
the experiment was sensitive to C. This is formulated mathematically 
by 




c - 






(111,44) 


where Ej^ < E 2 defines the extrapolated energy window, and 2 "^®) represents 

the differential energy spectrum for particles of mass A and charge Z, 

The correction factors computed as just described are listed in column A 

of table III. 3 together with the relevant energy windows at the detector 

and the top of the atmosphere. The energy window correction factors have 

also been computed on the basis of the assumption that all isotopes have 

12 

the same spectral shape as C. These are listed in column B of table 
III. 3. This has been done in order to facilitate comparisons with the 
data of other experimenters, who have generally made this assumption in 
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correcting their data (Webber et al. 1973, Preszler et al. 1975, Garcia- 
Vfunoz 1975a), However, we feel that the former approach is superior, 
and have used it exclusively in the interpretation of our data. We 
expect that the difference between the two approaches will' be most appar- 
ent when considering: 1) secondary components of the cosmic rays, whose 

spectra in interstellar space are expected to diverge from those of the 

12 

primary species (e.g. C) from which they are produced, and 2) isotopes 

which do not have A ~ 22 and are therefore affected differently by the 

12 

solar modulation process than is C. These observations are borne out 
by comparing the corre'ctions computed for ^^Be (a very important isotope 
which satisfies both conditions above) under the two sets of assumptions. 
The correction we have adopted has the effect of giving more ^^Be in 
interstellar space. 

3. Atmospheric Corrections . The production and destruction of 

the various isotopic species in the atmosphere above the detector has 

been corrected for using a model based on the semi-empirical formula of 

Silberberg et al. (1973) for the partial cross sections (to produce 

species j from an interaction where spacies i is incident on species k) . 

Equation III. 42 has been assumed for the total cross sections. The semi- 

empirical formula has been normalized to the results of Lindstrom et al. 

12 16 . , ‘ 

(1975a) for the particular cases of C and 0 primaries. The normal- 
ization factors adopted are listed in table III, 4. 

We- begin by temporarily treating the atmosphere as a slab of 
uniform thickness and known composition, the flux of species i surviving 


at depth X without undergoing any interaction is 


loCi) = 1,^0) e 


-Vi 


(III.45) 
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TABI£ III. 4 


Correction Factors Applied to Semiempirical Formula 
for Partial Cross Sections (Silberberg et al. 1973), as 
Derived From Data of Lindstrom et al. 1975a 


Priniarv 


Secondarv 


0.71 
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where 


. = Za -n I 

-Ur 


(III. 46) 


We can also write down a general formula for the flux of particles of 
species i at depth X which have undergone precisely k interactions: 


Q -r I •s'') 

x,,U)=J«'e ^ ^ 




(III. 47) 


where 




(III. 48) 


Equations III,45~48 may be integrated in an iterative fashion, beginning' 
with k = 1. The results, up to k = 3, are 


iLCi) =lto)e 


-■Si 2. 






e- ^ 


— 2^ J>(0) 




- ^i'Z. 

e ^ 


(III. 49) 




X 1- 




i 


^ ^,SL 

p <) p p 

^ i- t - — 

cv«.pi«i.-ftj)(Vf^P 
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We thus have, for the total flux of species i arriving at depth X after 
an arbitrary number of interactions 






(III. 50) 


In practice, the decrease monotonically with increasing k above 

-2 

some k'. For the present case (X ^ 5 ,g cm ), k' ^ 1, and terras beyond 

k = 3 are neglected. This is reasonable, since it is observed the J^2 

s; 0.02 J,* and J._ 0.0005 Jj . It should be noted that the J., (X) are 

^1.3 ^ ik 

all linear in the source abundances, J^(0), so that one may write the 
result III. 50 in the form 


J^{_o) CIII.51) 

Sp) 

Here the ik(X) are determined directly from equations 111,49,50. The 

s 

matrix D(X) may be regarded as the "propagator" matrix which propogates 
the cosmic-ray abundance vector through a slab of thickness X. Now, 

in reality the atmosphere is not a slab for our case. The experiment does 
not remain at a constant altitude, and particles are incident at various 
zenith angles. The actual distribution of atmospheric depths (effective 
slab thicknesses), f(X), is easily determined from the spark chamber data, 
and is shown in figure III. 13. This is readily incorporated into the 

O 

analysis by merely replacing D(X) by the atmospheric propagator matrix, 
■^D, defined by 

A C 

P = \ olX 4U) (111.52) 

where f(X) has been assumed to be normalized to unity. Then 




(III. 53) 



NUMBER OF EVENTS 


0.25 



Figure III. 13. Distribution of atmospheric slant depths. 
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The problem is now easily inverted to obtain the abundances at the top 
of the atmosphere in terms of those observed at the detector: 

This is obviously the formulation of interest here. 

A ‘1 

The general form of the matrix D is triangular. The diagonal 
elements are positive and greater than one to correct for attenuation in 
,the atmosphere, while the off-diagonal elements are negative to correct 
for contributions to the observed flux from spallation of heavier species. 
In order to compute the flux of species i at the top of the atmosphere, 

It is thus necessary to know the for all species as heavy 

or heavier (i.e, for 2 : Aj^) . Since we have only measured the isotopic 
abundances of the elements Be, B, C and N, it is therefore required to 
assume abundances for other species which might contribute. For this 
purpose, we have adopted the charge abundances measured by Ormes et al. 
(1975) using the same instrument in the CxS mode of analysis (i.e. Dl-3 
only) for E s 570 MeV/nuc. The relative abundances of the isotopes of 
each element were then taken from the predictions of Tsao et al. (1973) 
for the arriving cosmic rays. The various isotopes included in the 
calculation are indicated in table III. 5. 

In table III. 6, the isotopic abundances both at the detector 
and as extrapolated to the top of the atmosphere are given. The present 
data have been corrected for interactions in the detector (column A of 
table III. 3) and energy windows (column B of table III, 3) before perform- 
ing this calculation. The atmospheric correction factors listed in 
column D of table III. 3 are computed by merely taking ratios of the 
abundances listed in table III. 6. Once the Jj^(O) have been obtained. 
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TABl£ III. 6 

ZioCopic Abundances at Detector and Top of Atmosphere 



Detector 

Top of 


Detector 

Top of 


(~4.5rciii"2) 

Atmosphere 


(•4 ,5Rcm"^) 

Atmosphere 

’Be 

57-»4 

43+4 

*°Ca 

1042 

11+2 

"Be 

3D+4 

23+4 

♦iCa 

1+1 

1+1 

lOBe 

15+4 

icm 

*aCa 

4+1 

4+1 




*®Ca 

3+1 

3+1 

lOB 

93+10 

80+10 

**Ca 

3+1 

:i±i 

llB 

214+12 

197+12 







*«Se 

6+1 

7+2 

lac 

924+17 

932+17 




lie 

764B 

68^ 

48 Ti 

4+1 

4+1 




*7Ti 

4+1 

4+1 


U9+22 

111+23 

*®T1 

4+1 

4+1 

isu 

16C+17 

155+18 

49X1 

1+1 

l+L 


879+18 

924+19 


4+1 

3+1 

'’0 

23+3 

21+3 

80y 

2+1 

2+1 

18q 

18+3 

17+3 

Sly 

1+1 

1+1 

lay 

28+3 

26+3 

“®Cr 

2+1 

2+1 




siCr 

5+1 

4+2 

aoHe 

11446 

119+7 

s2Cr 

6+1 

6+2 

*iKe 

15+2 

13+2 

B3cr 

1+1 

1+1 

«Be 

26+3 

26+3 

8*Cr 

1+1 

1+1 

"®N« 

4244 

4244 

8a Hn 

6+1 

7+2 




8-*Ml 

4+1 

,4±1 

^Mg 

136+7 

147+8 

8«Mq 

4+1 

4+1 

»«Hg 

27+3 

27+3 




"««g 

28+3 

29+3 

64 fe 

5+1 

6+2 




8®Fe 

3+1 

3+1 

a’Al 

3644 

3844 

»"Fe 

58+4 

70+5 




“’Fe 

1+1 

2+1 

*esi 

11946 

131+7 




»Si 

ld+2 

10+2 

88Hi 

2+1 

3+i 

10 Si 

^2 

9+2 

eoui 

1+1 

1+1 

lip 

9+2 

9+2 




•as 

21+3 

23+3 




a»S 

4+1 

4+1 




a*S 

A 

5+2 

A^ndances and Errors 

for Be, B, C 

1«S 

1+1 

1+1 

and N 

from Maximum Likelihood Fit 




(Table 

II. 2), Other Charge Abundances 

» 8 C 1 

5+1 

5+2 

from Ormes et al. (1975) ; Isotopic 

** C 1 

3+1 

2+1 

Abundances from Silberberg and Tsao 




(1973) 

1 Errors Scaled 

from 27. for 

»«A 

6+1 

6+2 




»’A 

2+1 

1+1 




a*A 

5+1 

4+2 




«A 

1+1 

1+1 




»*K 

4+1 

4+1 




40k 

3+1 






3+1 

341 
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they way be propagated back through the atmosphere as is shown in figure 
III. 14, By using the abundances at 3 g cm instead of at the top of the 

atmosphere, we may. compute correction factors to an atmospheric depth of 

-2 

3 g cm , as given in column E of table III. 3, in order to facilitate 
comparisons with experiments performed at this depth. From table III. 6 
and figure 111,14, it is evident that the relative abundance of the 
important isotope is strongly affected by the atmosphere. The 

various percentage contributions to the ^^Be abundance observed at the 
detector are listed in table III. 7. 

In general, the percentage errors on the abundances are increased 
by the atmospheric correction. The major contribution to this increase 
often has nothing to do with uncertainties on the cross sections or 
abundances of other species. This may be seen by rewriting equation 
III. 54 in the form 


^ Ji ^ 



(III. 55) 


A -1 

Since the D are always negative (i.e. subtracting from the 

ik 

observed abundance what has been produced in atmospheric spallation 
reactions), we see that the percentage errors will be increased by the 
factor 

(III. 56) 

in performing the atmospheric correction. This factor approaches unity 
for species whose abundances at the detector are mostly primary (i.e. 
have not interacted in the atmosphere). For such cases. 




ABUNDANCE RELATIVE TO CARBON 


6 


iO ■ 



C ABSOLUTE 


0 10 
ATMOSPHERIC DEPTH (G CM “) 


Figure III. 14. Isotopic abundances relative to carbon 
as a function of depth in a slab model 
atmosphere . 
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primary 

Isotope 


10 

10 , 

11 . 


Be. 



14|j 


0 

F 


Ne 

Na 

Mg 

A1 

Si 

P 

S 

Cl 

A 

K 


Ca 

Sc 

Ti 

V 


Cr 


■Mn 


Fe 

Co 


Ni 


TABI£ III. 7 

CONTRIBUTIONS TO OBSERVED 

Number of ^^Be per 10,000 Observed Originating 
in Atmospheric Interactions of Heavier Species 

6413 

0 

1118 

741 

84 

60 

175. 

594 

26 

132 

39 

160 

34 

107 

8 

28 

8 

12 

11 

23 

10 

22 

9 

23 

26 

131 

0 

6 
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XU) » J.U) <I«.57> 

Errors on the abundances of all species contributing to an observed 

abundance have been considered. In cases where the abundances at the 

detector have been inferred from the charge abundances measured by Ormes 

i 

et al. (1975), the errors are taken to scale as with a 2% uncertainty 
for ^^0 as the normalization. The uncertainties on the abundances both 
at the detector and as extrapolated to the top of the atmosphere are 
indicated in table III. 3 and also table III. 6, 

4, Demodulation Corrections . Before any meaningful comparison 
of the results with model predictions can be made, the data must be cor- 
rected for the effects of solar modulation. Because of the large uncer- 
tainty involved in making this • correction, it is not attempted in comparing 
with other experiments. It is only used for comparing -the data with model 
predictions as discussed in Chapter IV. The procedure by which this 
correction is computed is completely analogous to that used for the energy 
correction discussed above. The model of Garcia-Munoz et al. (1975b), 
which was used to predict the local interplanetary spectra required for 
computing the energy corrections, also predicts spectral shapes in local 
interstellar space. Provided one can map the energy windows in which 
the experiment is sensitive to- local interstellar space, relation III. 44 
may then be applied again with the substitution of interstellar spectra 
for the interplanetary spectra. The demodulation correction would then 
simply consist of the ratio of the correction to local interstellar space 
to the correction to local interplanetary space (i.e. the previous energy 
correction) . 
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In order to extrapolate the energy windows to local interstellar 
space, it is required to understand the process of solar modulation. 

This process becomes increasingly important as the energy of a particle 
is decreased below ~ 1 GeV/nuc, and is thought to result in the virtual 
exclusion of particles below ~ 200 MeV/nuc. The process is not completely 
understood at the present time. Almost all current models, however, are 
derived from one first introduced by Parker (1958). In this approach, 
the modulation process is modeled in terms of a one dimensional diffusion 
of the cosmic rays into the solar system along magnetic field lines being 
convected outward from the sun by the solar wind. The quantitative form- 
ulation of the model is given by (see Parker 1958, 1963; Gleeson et al. 
1968; Fisk et al. 1969 for further details) 


\ 


(III. 58) 


where 




1 jEi- wc 

0 E -t Tnc 


(III. 59) 


i« the number density of cosmic rays of species i, V - 400 km sec 
is the velocity of the solar wind, r is the heliocentric distance, K is 
the parallel diffusion coefficient and E is the kinetic energy. Here the 
first term describes the convection of the cosmic rays in the solar wind, 
the second describes diffusive effects due to scattering of the cosmic 


rays from inhomogeneities encountered along the magnetic field lines, and 


the last term accounts for adiabatic energy losses in the expanding solar 
wind. This model deals only with motion confined entirely to the equator- 
ial plane of the solar system, and thus ignores parameters which are 
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thought to be important, but have not yet been investigated experimentally 
(Cecchini et al. 1975, Fisk 1975, Moraal et al. 1975). Equation III. 58 
has been solved numerically by Fisk (1971). However, if the streaming 
of cosmic rays 





1 - 




ilL 




(III. 60) 


vanishes, then it is possible to obtain an analytic solution to 111,58. 
This is known as the force field solution, and is due to Gleeson and 
Axford (1968). The result may be formulated as 


= jfCPb, n.) 

where 


(III. 61) 


(III. 62) 


relates the differential rigidity and energy spectra, and 

('^ Kf 2 ■£ (III. 63) 

Jr^ 3\<y- - \ £ -f. 7 t)c^ 

(with K = assumed) relates r^ and r^. Thus the differential 

rigidity spectra at r^ and r^^ are simply related to each other. Garcia- 
Munoz et al, (1975b) at the University of Chicago have used the numer- 
ical technique of Fisk (1971) to obtain forms of the diffusion coef- 
ficient, K, and interstellar differential energy spectra, J^, which are 
consistent with the spectral shapes observed near the earth. The result 
is 




for 


X cX ^ 4oo HeV) 
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-ftfr )r> 50 /^U 
4^ Y-4 ^0/lU 

( 111 . 64 ) 

£ > 0.4 

( 111 . 65 ) 


at the cosmic-ray source. Using equation 111.64 for the diffusion co- 
efficient, equation III. 63 may be integrated to give the mean energy 

loss of particles penetrating from outside the solar cavity (r -* <o) to 

b 

the orbit of the earth (r ~ I AU) : 

a 


Efc-Ea, ^ 2 MtV 


(III. 66) 


We emphasize that equation III. 66 gives only the average energy loss 
experienced by a particle of given Z, A and g in penetrating the solar 
cavity, and cannot be meaningfully applied to individual particles (e.g. 
see Goldstein et al. 1970). 

The demodulation correction factors computed as described above 
(i.e. using the force field approximation to map the energy windows to 
the local interstellar space and assuming the interstellar spectral shapes 
calculated by Garcia-Munoz et al. 1975b) are listed in column F of 
table III. 3, along with the extrapolated interstellar energy windows. 

Due to the uncertainty introduced by the process of solar modulation, 
the local interstellar spectra are not well known at low energies. The 
best one can say with any certainty, is that they probably lie somewhere 
between pure power laws in rigidity and pure power laws in total relativ- 
istic energy per nucleon (Ramadurai et al. 1972). It is thus desirable 
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to see the latitude which might be introduced by the assumption of var- 
ious possible spectral shapes on the correction factors applied to the 
observed isotopic abundances. This has been done in table III. 3 for the 
cases of interstellar spectra which are pure power laws in rigidity 
(total relativistic energy per nucleon, or Lorenz factor) with spectral 
index -2.60 — column 1(G) — and spectral index -2.75 — column J(H). These 
can be compared to column F, which gives the correction factor we have 
actually used in correcting the data. This correction is based on more 
realistic forms of the interstellar spectra, and in particular takes 
account of differences in the spectra of primary and secondary components 
of the cosmic rays. Evidently, the various demodulation corrections all 
affect the isotopic abundances corresponding to a given element in 
roughly the same way. The major difference is in the pronounced charge 
dependence of the correction factors based on the assumption of rigidity 
power laws. 

F. Discussion of Results. 

The results of the present experiment are essentially contained 

in tables III. 2 and III. 3. Although the resolution is not quite good 

enough to see a separate peak for ^^Be, it is comparable to what has 

10 9 

been reported by other experimenters. Since Be and Be are not clearly 
separated in the data, a fitting procedure has been relied upon to deter- 
mine the relative abundances of Be isotopes. Furthermore, the uncertainty 
in the abundance is increased from 25% to 40% by the presence of the 

atmosphere above the detector. The factor limiting the resolution is 
probably the detector normalization. This is closely related to the 
necessity of using semirelativistic particles to determine detector 
normalizations, and also differences in the scintillator response functions 



113 


from detector to detector. A detailed accelerator calibration of the 
instrument would thus be of immediate value in improving the resolution 
obtained. 

G. Comparisons with Other Experiments. 

In the following section, the results of the present experiment 
are compared with those of other experimenters. Experiments using the 
geomagnetic cutoff technique (Juliusson 1975, Dwyer et al, 1975) have 
been excluded from the comparison for the following reasons: 1) the 

geomagnetic cutoff experiments are generally performed at substantially 
higher energies (2 0.7 GeV/nuc) than is the present experiment 0.2 
GeV/nuc); 2) the cutoff experiments so far have not yielded more than 
mean masses for the various elements although it is in principle pos- 
sible to derive the actual isotopic composition (Peters 1974); and 3) 
the assumptions which go into the cutoff technique (especially with re- 
spect to energy spectra) are at best approximations, the reliability of 
which make the interpretation of the results somewhat uncertain. Insofar 
as a comparison ^ possible, however, there are no outstanding dis- 
crepancies between the present results and those of experiments using 
the geomagnetic rigidity cutoff technique. 

In table III. 8, we have compared the present data with the 
results of other experiments. Included are the pioneering results of 

Webber et al. (1973), which are balloon observations carried out at an 

_2 

atmospheric depth of ~ 3 g cm and about the same energy as the present 
experiment. This group has also recently reported results of an improved 
experiment (Freszler et al. 1975), which is also included in the table. 
Finally, the IMP 7 and IMP 8 results of Garcia -Munoz et al. (1975a) are 
also included. These are satellite observations taken outside the 
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Exparlatnc 
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MlAK-2 

CSPC 

CHIC-l 

CHIC-2 

Correct Iona* 

BCE 

... 

— 

BCD 

— 

— 

Energy (MeV/m*c) 

--250 

-2C0 

-200 

-250 

-100 

-100 

Ac»»pherlc 







Dcpth(Rcc'2) 

(3.0) 

2.9 

3.1 

(0) 

0 

0 

Abundance Ratio: 
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0.57 -*0,05 
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0.60 

0.63 40.07 
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0.93 t0.03 
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0.07 Vo, 01 



0.07 40.01 



C/C 

l.C00^0.026 


1.000 

1 .0402.0.036 

1.00040.014 

1,00040.017 

»*H/H 
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0.40 

0.43 *0.09 




0.58 40. C7 


0.60 

0.67 4^0.08 
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0.27540.017 


0.271 

0.7902.0.026 
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0.22840.008 
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atmosphere and magnetosphere of the earth, at energies roughly half those 

1 

at which the present data are taken. The results of these various experi- 
ments are referred to by NHAM-1, NHAM-2, CHIC-1 and CHIC-2, respectively. 
In each case, our results have been corrected in such a manner as to be 
most consistent with these observations at the respective detectors. 


In particular, the present data have been corrected to an atmospheric 

depth of 3 g cm (column E of table III. 3) for comparisons with NHAM-1 

and NHAM-2, and to the top of the atmosphere (column D of table III, 3) 

for comparisons with CHIC-1 and CHIC-2. In all cases, the data have 

been corrected for interactions in the detector (column C of table III. 3) 

and for energy windows assuming that the energy spectra of all isotopes 

12 

have the same shape as that of C (column B of table III. 3). The agree- 
ment among the various experiments considered in table III, 7 is generally 
quite good. The only significant discrepancy among the relative isotopic 
abundances reported concerns the isotope ^^Be, All of the results 


10 

reported are consistent with Be s: 0,6, except for those of CHIC-2, 

10 

which ate most consistent with the complete absence o'f Be. The low 
10 

^^Be abundance observed by CHIC-2 is compensated by the largest Be 

Be 

ratio of any of the experiments reported. All the remaining discrepancies 
concern the relative charge abundances. These are as follows: 1) CHIC-1 

and CHIC-2 both obtain a smaller value of the ratio Be/C than does either 
NHAM-2 or the present experiment; 2) the present experiment observes 
the largest B/C ratio of any of the results reported; and 3) CHIC-1 
and CHIC -2 both observe a smaller N/C ratio than does either NHAM-2 or 
the present experiment. Finally, we note that the present results agree 
reasonably well with other results on the charge composition obtained at 
similar energies (e.g, see Shapiro et al. 1973). There is also a reason- 
able agreement with the results of Ormes et al. (1975), which were 
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obtained from the same experiment (in the C x S mode of .analysis) at 
somewhat higher energies 570 MeV/nuc) and with superior statistical 
weight. 


Let us now consider the degree of consistency of the results 

presented in table III. 8. In particular, we will concentrate on the 
10 

ratio which represents the only significant discrepancy among the 

Be 

relative isotopic abundances compared. 

We first note that all of the results except CHIC-2 would be con 
10 

sistent with the value 2 : 0.06 at the top of the atmosphere. Of 

the other results, the largest deviation from this value is by the present 

experiment, which is one standard deviation high. We recall from the 

10 


discussion in section C.5 that our observed Be abundance may be reduced 

by about one standard deviation if we assume that some of the observed 

10 9 

Be is due to an asymmetry in the Be mass distribution. Although such 

10 

an effect would considerably reduce the scatter in — among experiments 
other than CHIC-2, it cannot by itself make the present result consistent 
with the complete absence of in the cosmic rays. The above consider- 

ation represented the only reasonable possibility for error in our observed 

10 „ 

value of — Other possibilities which were considered not probable 
Be 

included a breakdown in the mass scaling of the variance selection 

threshold and the standard response curves, and an error in the atmospheric 

correction to the data due to the use of erroneous cross sections. Thus 

it is quite difficult to see how the present experiment could be made 
10 ^ 

consistent with be ^ 0 . In particular, only an improbably large error 

Be 

in the atmospheric correction could at once make the present results 
and those of NHAM-1 and NHAM-2 consistent with the absence of ^^Be. 


However, this seems to be very unlikely. 
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We next observe that, while the Chicago group observes the smallest 
10 ^ 

value of 5^ at the lowest energies 100 MeV/nuc), the present experi- 

Be 

ment observes the largest value at the highest energies 250 MeV/nuc). 

Thus there exists the possibility of a real energy dependence of the 
10 

ratio __E^. One possibility for understanding such an energy dependence 
Be 

can be seen from the assumption that the mean pathlength traversed by 

cosmic rays, is a constant independent of energy, at least below 

1 Gev/nuc. In this case, equation I.l predicts that the mean confinement 

time of cosmic rays in the galaxy, Tg, will vary as 1/p. A second effect 

which may contribute concerns the energy dependences of the various 

cross sections for the production of from interactions of heavier 

species with the interstellar medium. Both of these points have been 

taken into account in the propagation calculations of Garcia-Munoz et al. 

(1975b) which predict a decrease in the ratio of about 207o in going 

from 250 MeV/nuc to 100 MeV/nuc. This -of course assumes a mean energy 

loss of ~ 250 MeV/nuc per particle in penetrating the solar cavity. If 

lOg 

we disregard solar modulation, then the ratio — g— changes by 50% between 

250 and 100 MeV/nuc. Clearly even this is inadequate to explain the 

observations. Thus if there exists a real energy dependence of the ratio 

lOv, 

■ j> , then there appears to be no ready explanation for its magnitude. 

Before concluding, we wish to make a few general remarks with 
respect to the various experiments. Both the present experiment and those 
of the University of New Hampshire group are balloon borne experiments 
performed under a residual layer of the earth's atmosphere. The IMP 7 
and IMP 8 experiments of the Chicago group, on the other hand, are per- 
formed well outside the atmosphere and thus avoid any effects of the 
atmosphere on the observed data. However, we have seen that the cross 
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sections used to make the atmospheric corrections to the present data 

are not likely to be sufficiently in error to have a significant effect. 

Thus the principle disadvantage imposed by the presence of the atmosphere 

10 

is to increase the uncertainty of the observed Be abundance from 25% 
to 40% after extrapolation to the top of the atmosphere. On the other 
hand, the present experiment and those of the University of New Hampshire 
group are performed at significantly higher energies than are those of 
the University of Chicago. Although this does not directly affect the 
measurement, it does significantly affect its interpretation. This 
follows since a highly uncertain correction for the effect of solar 
modulation is required before a meaningful comparison with theory is 
possible, and the uncertainty of this correction increases with decreasing 
energy. Finally, we feel that the means by which background is removed 
in the present experiment is superior to that used by either of the 
other groups » This derives essentially from the fact that only the pre- 
sent experiment is capable of "continuously" monitoring the identity of a 
particle over its entire trajectory in the experiment. 



CHAPTER IV 


DISCUSSION OF RESULTS 

Once the experimental results have been derived and corrected to 
yield isotopic abundances in the local interstellar medium, comparisons 
may be made with the predictions of the various models for the propagation 
and origin of cosmic rays. Although some models have been proposed in 
which the cosmic rays fill some extended region of extragalactic space 
(e.g. see Setti et al. 1971, Brecher et al. 1971 and Sitte 1972), most 
current models propose the confinement of the cosmic rays either to the 
galactic disk or some extended region (known as the galactic halo) sur- 
rounding it. The actual existence of the galactic halo is highly con- 
troversial, and it is an objective of the present thesis to definitively 
answer the question of whether it is required in order to confine the 
cosmic rays . 

A. The Diffusion Model 

Halo Model . Most currently popular models for the propagation 
of cosmic rays in the galaxy are based upon the diffusion model first 
introduced by Ginzburg (Ginzburg et al. 1964). The essential observation 
is that although cosmic rays of all but the highest energies (i.e. E ^ 

10^^ eV/nuc) are constrained to closely follow magnetic field lines 
through the interstellar medium, these field lines are sufficiently 
randomly oriented in space that even a very small drift (or scattering) 
across the field lines (as would be caused by irregularities encountered 
in the magnetic field) essentially randomizes the trajectory of a particle 
through space. Ginzburg postulated that the cosmic rays originated in the 
galactic disk, probably as a result of supernova events, then diffused 
freely throughout a quasi-spherical region surrounding the disk known as 
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the galactic "halo". Although the magnetic field intensity and gas 

-8 -2 "3 

density postulated to exist in the halo (~ 10 gauss and ~ 10 cm , 

respectively) are two orders of magnitude less than those characteristic 

-6 ••3 

of the disk 10~ gauss and 1 cm ), particles were imagined to both 
escape from the disk into the halo and re-enter the disk freely. Ginzburg 


.29 


-1 


found that a diffusion coefficient given by Kjj s; 10 cm sec to be 

appropriate for such a model. This results in an effective scattering 

3Kx, 




3pc, and a mean escape time for cosmic 


mean free path given by Xy 2 ? 

^ 8 

rays from the halo given by Te s? 2 x 10 yr (where r„ 12 kpc) is 


2K 


■H 




the halo radius. This gives a reasonable value for the mean pathlength 

_3 

traversed by cosmic rays (Xg s; 5 g cm ) as deduced from the abundances 
of the elements Li, Be and B (presumed to be absent in the source) 
relative to their progenitors, such as C and 0. In addition, the scat- 
tering mean free .path obtained is of the same order as the inhomogeneities 
observed in the galactic magnetic fields (Allan 1972). Finally, the 
cosmic rays are scattered sufficiently before observation (~ 10^ times) 
to lose all memory of their points of origin (i.e. source locations), 
which is important in order to reproduce the high level of isotropy 
observed in the cosmic rays. A further advantage of such a model is 
that it considerably reduces (by a factor ~ 10^) the time required for 
particles to escape from the system, and thus considerably reduces the 
energy requirement on the cosmic-ray sources (i.e. the rate at which 
they must produce cosmic rays in order to replace those escaping and, 

thus, maintain the cosmic-ray density at a constant level). In quanti- 

-3 

tative terms, for the observed cosmic-ray energy density w s; 1 eV cm , 

CR 

40 -1 

the power required of the cosmic-ray sources is s; 5 x 10 erg sec 

68 3 

for a halo volume of ~ 2 x 10 cm . For comparison the power available 
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from supernova explosions is Pgjr ^ 10^^ erg sec“^, assuming a supernova 
r-ate o£ < 1/50 years and < 10 erg/supernova in cosmic rays produced. 

The power required of the cosmic-ray sources in the case where cosmic 
rays freely escape from the galaxy is ~ 10 times that required in the 
Ginzburg diffusion model, as may be seen by comparing Tg and Thus 

it is clear that supemovae can be the primary source of energy for cosmic 
rays in the galaxy only if there is some mechanism available to restrict 
their escape, as (for example) diffusion in the Ginzburg model. This is 
important since supernovae are among the most powerful local sources of 
energy we can imagine. 

It is apparent from the foregoing discussion that the model of 

Ginzburg, which postulates diffusion of the cosmic rays in a galactic halo 

meets the main requirement of a cosmic -ray propagation model. It is 

possible to reproduce the characteristic mean pathlength traversed by 

_2 

cosmic rays (\g ~ 5 g cm ) and also their high degree of isotropy 

(6 ^ 10”^ for E > 100 GeV/nuc; Allan 1972). At the same time, the 

cosmic rays are confined In the galaxy sufficiently long that an unreason 

able burden is not placed on the sources. This is all accomplished with 

a diffusion coefficient ~ 10 cm sec” in a halo of radius r^ ~ 

-2 q 

12 kpc with gas density njj ~ 10 cm and magnetic field intensity Bjj ~ 
10“° gauss, all of which is quite reasonable. Halos have been observed 
around various other galaxies (not by any means all, however) by means of 
the synchrotron radio emission of the cosmic-ray electrons trapped in 
them (Baldwin 1954). The existence of a similar halo around our own 
galaxy, however, is by no means firmly established and is, in fact, 
strongly questioned by radio astronomers (Burke 1967). 

2. Disk Model. The disk model, a variant on the halo diffusion 
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model o£ Ginzburg, would confine the cosmic rays predominantly to the 

galactic disk (Owens 1975; Dickinson et al. 1975), and thus circumvent the 

problem imposed by the failure to observe a galactic radio halo around our 

own galaxy. We note that this model does not preclude the existence of a 

cosmic-ray halo, but merely postulates that the re-entry of cosmic rays 

into the disk from the halo is negligible. Both the gas density. and the 

“3 -6 

magnetic field intensity in the disk (n^ ~ 1 cm and ~ 10“ gauss, 

respectively) are about two orders of magnitude greater than the values 
characteristic of the halo. The former consideration requires an escape 
time two orders of magnitude smaller than that in the halo model (i.e. 

S: 2 X 10^ yr) if the observed mean pathlength traversed -by the cosmic 
rays s: <nt>njjcTg ^ 5 g cm ‘^) is to be maintained. Since most particles 
will reach the top or bottom face of the disk before reaching its circum- 
ference, the half thickness of the disk (Zj^ s: 0,25 kpc) replaces the halo 

radius in relating the escape time to the diffusion coefficient: 

3z 2 23 2 — 1 

. P . We thus deduce a diffusion coefficient, s; 10 cm sec , which 

2^D 

is an order of magnitude smaller than that appropriate to the halo model. 

This corresponds to a correspondingly smaller effective scattering mean 
3K 

free path, \ ^ — £ = 0.3 pc. In this model, the power required of the 

XJ o 

cosmic-ray sources, Pqrs s: 1.5 x lo"^^ erg sec ^ is only three times that 
required in the halo model, and still within reason. 

The halo and disk diffusion models for the propagation of cosmic 
rays which have been discussed above characterize the two extremes of a 
continuum of models in which the halo is gradually flattened, and eventually 
merges into the disk. All of these models can be made consistent with 
observations. Perhaps the most sensitive test of which model is correct 
would be the direct determination of the mean matter density, n, in the 
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region where the propagation takes place, or alternatively the mean escape 
time, Tg, to which it is inversely proportional (by virtue of the known 
leakage pathlength, 2= <nt> nc Tg). The radioactive secondary isotope, 
^^Be, with a mean decay time, = 2,2 x 10^ yr, provides an ideal tool 
for doing this. 

Other variations on the basic diffusion model have also been 
proposed. These include the assumption of boundary layers (i.e, for the 
disk, halo, etc.) which are either partially or totally reflective. Such 
boundaries may, for instance, be set up by self-generated magnetohydro- 
dynamic waves. One may also replace the continuous distribution of 
sources in space and time by a more realistic distribution of discrete 
sources, which is especially important in the case where one of the dis- 
crete sources happens to be quite nearby. Diffusion in only one dimension 
(i.e. along magnetic field lines) may also be considered, or alternatively 
"compound diffusion" in which one considers separate diffusion processes 
for the random walk of magnetic field lines in space and the motion of 
particles along the field lines. As a final example, one may consider 
multiple confinement regions, each of which is characterized by a dif- 
ferent magnetic field intensity and gas density. In particular, supernova 
envelopes with characteristically high gas densities and magnetic field 
intensities may well surround the cosmic-ray sources. The cosmic rays 
would then be required to escape from these before entering the galactic 
disk, 

B. The leaky Box Approximation 

A rather general formulation of cosmic-ray transport in the dif- 


fusion model is given by 
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(iv.l) 


where J^Cpj'Tjt) represents the differential flux of cosmic rays of 
species i and momentum "p at position r and time t. Qj_Cp,?,t) represents 
the source function for species i, and includes production and loss in 
nuclear interactions and radioactive decay, and is the diffusion 
tensor for species i. If the cosmic-ray flux is isotropic, we may then 
replace ^ by 7 everywhere in. equation IV. 1, noting that 

(IV.2) 


relates the differential momentum and energy spectra. It is usual to 
replace the diffusion tensor by a scalar since there are no distinguish- 
able directions in the problem. Furthermore, the diffusion coefficient is 
often taken to be independent of 7 , r and t. To the best of our knowledge, 

the cosmic-ray intensity, has been constant over at least the last 
9 

~ 10 yr (Van Loon 1973), so that it is reasonable to formulate a steady 
state model by neglecting the first term of equation IV. 1. Consider next 
the second term. For 7 ^ 1.5, we may put 


^ /V. 

- 


f/3c[6-M^V fcjri) 


ft 



^ Li* JO ^)fc /) 


(IV. 3) 


Assuming energy spectra which are power laws in total energy per nucleon 
of spectral index ~ -2.7, we have 


ar 



fcA7^ 


r 


(IV. 4) 
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For comparison, one term which contributes to is the loss to nuclear 

‘interactions. Taking the total inelastic cross section to be ^tot — 

2/3 

50 A mb, this term is 

f C. (IV. 5) 

We thus see that the second term of equation IV. 1 is negligible for 7 > 
1.5, at least at the 10% level. After applying all of the approximations 
discussed above, equation IV. 1 is reduced to 









(IV. 6) 


which is a commonly used formulation of the diffusion model. We have 
already seen in section II. A that the mean time for cosmic rays to escape 
from the system in the halo (disk) model is 2 x 10 yr (~ 2 x 10 ° yr) , 
which corresponds to a total distance traveled of ~ 60 Mpc (~ 0.6 Mpc), 
This is ~ 5000 (2000) times the distance to the nearest boundary, and 
suggests a picture in which the cosmic rays approach the boundaries often 
in their lifetimes with a low probability of escape each time. A 
reasonable approximation in this case is to neglect the positional depend- 
ences and integrate equation IV . 6 over all space interior to the 
boundaries : 


f K. av.7i 

Here the surface integral. is interpreted as the rate at which particles 
cross the boundary surface, and is approximated by the last equality. 
Here Jj^(7) Is the mean differential cosmic-ray intensity Inside the 
boundary and is the mean time it takes a typical cosmic-ray particle 
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to leave the bounded volume. We may now write the result 


Q..(r) = (IV.8) 

This is known as the leaky box model of cosmic-ray propagation. It may 

be interpreted in terms of a leaky box (i.e. the boundary surface) in 

which the cosmic rays are confined. The probability that a given cosmic 

dt 

ray will escape from the box (i.e. leak out) in time dt is merely — , 

Te 

and this loss is exactly compensated by the production of Q^( 7 )dt, new 
cosmic rays, either directly by sources or in the spallation (or radio- 
active decay) of heavier species. The leakage time, may be related 
to the diffusion coefficient, K, according to 



VcW ^ * 

) 2f.\o jr 


ynocfjJL 

(IV. 9) 
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TnoJ-ell. 



where is the net distance a particle must propagate from its 

source (on the average) in order to cross the boundary. We have assumed 
that the diffusion coefficient K (and hence also t^) is independent of 
species. We proceed by expanding the source term, Q^(7), in equation 
IV.8 to show explicitly the contributions of production and' loss in 
nuclear spallation reactions and radioactive decays; 





t 


(IV. 10) 


Here denotes that part of the source term which is directly due to 
cosmic -ray sources, and specifically does not include nuclear spallation 
or radioactive decay. < denotes an average over the target atoms in the 
confinement volume (V), n is the number density of such targets and’ 
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J * -1 

-* i) « is required. The divergence of the term (i ->•* ) 

*J 

in equation IV. 10 then requires either that J,(7) -► 0 with *+• *) 

or that become large in order to balance the divergence of this term. 
The latter alternative is clearly unphysical except as an extreme assump- 
tion. Consequently, J^(7) -»• 0 and species i does not contribute signif- 
icantly to any other species except the one to which it decays, species j. 
Thus species i may be absorbed into species j by the substitutions 


and 






jjyn 


q' 
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S.^- + 0-x. 


(IV. 13) 


(IV. 14) 


Radioactive isotopes which do not clearly fall into either of the two 

categories described above (i.e. long or short mean lifetime against 

decay) are most rigorously retained as they appear in the calculation. 

However, it often happens that such isotopes contribute negligibly to the 

k 

production of other species (through the partial cross sections, ) 

because of their relatively low abundances. In such cases, T^(j i) 

-»■ CO may still be an adequate approximation for the calculation of the 
source abundances, Q^(7), from equation IV. 11 as a function of the leak- 
age pathlength, Xg. The calculation may then be inverted using the actual 
value of ''’^(j -* i) iu order to determine the gas density, n, in the 
propagation region, or equivalently the mean escape time, . This will, 
in fact, be the approach taken below for the case of ^^Be. It should be 
noted that the approximations of long or short lifetimes of radioactive 
isotopes against decay as described above are not required in order to 
apply equation IV. 11 to the data. In particular, they do not serve to 
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eliminate from the problem any unknowns which would otherwise prevent its 

solution. The main advantage of the approximations lies in the resulting 

simplifications, particularly in the number of species which must be 

LB 

considered in constructing the leaky box propagator matrix, D. 

1. Path length . The application of the leaky box model to the 

interpretation of experimental data is now a simple matter. One has only 

LB 

to calculate the leaky box propagator matrix, D, for assumed values of 
the leakage lifetime, T^, and the density of matter in the confinement 
region, n. Given the observed abundances of the arriving cosmic rays, 

the source abundances, q|( 7), may then be calculated. The values of 
and n are then adjusted so that the abundances of all components assumed 
absent In the source actually do vanish. Examples of such isotopes would 
be those of Li, Be and B (which cannot survive in the hot environment of 
a star), and perhaps also isotopes such as C and which are not 

easily produced in the chains of nuclear reactions thought to be relevant 
to the cosmic-ray nucleosynthesis. 

Various authors have applied the leaky box model, essentially as 
formulated above, to the interpretation of the observed charge composition 
of the cosmic rays (e.g. see Shapiro et al. 1973). In such cases, it is 
also required to make assumptions regarding the isotopic composition of 
the various elements in the cosmic rays, either at the source or as 
observed. This is usually done using the observed solar system abund- 
ances (Cameron 1973) as a standard. The results are insensitive to the 
survival or decay of the few radioactive isotopes for which it is question- 
able, so that only the leakage pathlength, is determined with any 

9 -2 

precision. The results range from > ~ 3 g cm to \ ~ 7 g cm at 

e e 

moderate energies of a few hundred MeV/nuc to a few GeV/nuc. Results at 
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f^(j i) denotes the mean lifetime against radioactive decays taking 
s.pecies j to species i (* denotes "anything”)* As formulated in equation 
IV.lOj the leaky box models have two characteristic parameters which may 
be adjusted in order to achieve the optimum agreement with observations: 
the leakage time, and the number density of target atoms in the "box", 

n. The relative abundances of the various target atoms are assumed to be 
sufficiently well known (i.e. ~ 907« H and 10% He by number) as to intro- 
duce no extra degrees of freedom, ‘ Fpr stable isotopes not subject to 
radioactive decay, = pRcT^ = <m>^c(nTe) is the only free parameter 
(<n^ is the mean mass of the target atoms). We note that equation IV. 10 
is a matrix equation for the source abundances, Q^(7), in terms of the 
observed cosmic-ray abundances, J,(7), i.e. 



defines the leaky box propagator matrix. 

The application of equations IV. 10 through IV. 12 to the interpre- 
tation of experimental data is often facilitated by various assumptions 
regarding the stability of the radioactive isotopes. The terms involving 
'Td (j ->■ i) may be altogether neglected in equations IV. 10 through IV. 12 
if the lifetime against radioactive decay is sufficiently long. The 

condition for this approximation to hold is Tjj(j -> i) » Tg if S 

Ta 

otherwise, (j i) » ~ Tg is required. If the radioactive 

J • 

1 , . ' 

decay is very rapid, on the other hand, one may put T,j(j -♦ i) 0. The 

condition for this to hold is (j i) « if J,(7) ^ J. (7) 5 otherwise 

J ^ 
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higher energies seem to indicate a systematic decrease in X as one goes 

e 

to progressively higher energies (Lezniak et al. 1975, Juliusson et al. 
1975). 

The results of a similar calculation, as applied to the present 

results, are shown in figure IV. 1. Here the cosmic-ray source abundances 

are computed as a function of the assumed value of the leakage pathlength, 

X , on the assumption that ^%e is stable. The cross sections used in 
e 

this calculation are based on the same formulae as were used in computing 
the atmospheric corrections (section III.E.3), The cross sections have 
been computed based on an assumed propagation energy E = 600 MeV/nuc, 
and assuming the composition of the interstellar medium to be 10% helium 
and 907. hydrogen by number. The assumptions regarding the stability of 
the various isotopes considered in the calculation are given in table 
IV. 1. For each isotope, the range of source abundances consistent with 
zero has been determined based on the statistical errors assigned to the 
observed isotopic abundances based on the maximum likelihood method. 

This is indicated by the darkened portion of the curve for each isotope. 
The dashed extensions to these error limits also take into account 
estimated uncertainties in the cross sections used in the calculation, 
estimated to be ~ 10% for total cross sections and ~ 15% for partial 
cross sections. This has been carried through the corrections for inter- 
actions in the detector and the atmosphere as well as the leaky box 
propagation. The best value of for each isotope (assuming it is 
entirely secondary in origin) is that at which its source abundance is 
exactly zero. The error limits on X^ are then estimated from the range 
of values of for which the calculated source abundance is consistent 
with zero, i.e. that range over which the curves in figure IV, 1 are 
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Figure IV. 1. 


Isotopic abundances at the cosmic-ray source 
relative to carbon as a function of mean 
leakage pathlength, Xg, assumed. Be assumed 
stable. Heavy portions of curves indicate 
values consistent with zero if only uncertain- 
ties on isotopic abundances are considered. 
Dashed e.xtensions include effects of cross 
sect'fon uncertainties assumed to be 10% for 
total and 157o for partial cross sections. 
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darkened. The resulting estimates, of the leakage pathlength for each 
isotope measured are tabulated in table IV. 2, and illustrated in figure 
IV. 2. These results are compatible with those of other authors, and 
apparently indicate a real discrepancy between the values of the leakage 
pathlength obtained based on the Be isotopes (i.e. ^ 3.5 g ctn”^) and 

those obtained using other isotopes (i.e. ^ 6 g cm ). We also note 

that the only isotopes which significantly affected by the survival or 
decay of ^®Be are ^^Be itself and, to a much lesser extent, Thus 

the leakage length obtained for ^^Be is increased if part of this isotope 
has decayed, while that for Be is somewhat decreased. Finally, the 
contribution to the total error in the pathlength determination from 
uncertainties in the cross sections used is observed to be generally less 
than or about equal to that from the uncertainties in the measured 


isotopic abundances. 


2. Lifetime . In order to estimate the density of matter, n, in 
the cosmic-ray confinement volume consider equation IV. 10 for the specific 
case of ^*^Be (denoted by putting i = 10): 


<m> : 

■= ^ 


(IV. 15) 


i <m> ''r 

Here we have substituted = <nc> ncR/Tg, taken Q^q(/) =0, and noted 
explicitly that no other isotopes decay into ^^Be, Observing that the 
right hand side is independent of the survival or decay of ^^Be, we may 
write equation IV. 15 again for the case where does not decay: 


li. r 


(IV. 16) 



TABLE IV. 2 


Leakage Pathlengths Derived from Various Isotopes 

(g cm“^ units) 

Lower Limit Upper Limit 

Best 


Isotope 

B 

A 

Value 

A 

B 

^Be 

3.28 

3.49 

4.01 

4.53 

4.88 

^Be 

2.46 

2.58 

3.24 

4.04 

4.27 

l°Be* 

0.89 

0.91 

1.60 

2.40 

2.46 


4.09 

4.37 

5.50 

6.48 

7.07 

11b 

4.93 

5.54 

6.33 

6.88 

8.02 

12„ 







13c 

4.70 

5.10 

6.53 

7.89 

8.87 

l^N 

7.61 

8.22 

— 

— 

— 

15n 

8.12 

9.53 





_ . _ _ 


A: No Cross Section Errors Included 

B; Cross Section Errors Included in Interaction Correction, 
Atmospheric Correction and Leaky Box Calculation 

*10 


Be assumed Stable 
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ISOTOPE 

Figure IV. 2. Mean leakage pathlength, as derived from con- 
dition that source abundances of various isotopes 
vanish. Dashed extensions to error limits include 
contributions of cross section uncertainties assumed 
to be 10^ for total and 15% for partial cross sections. 
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We Tnay now eliminate the right hand sides of equations IV. 15 and IV. 16, 
and solve for the density, n. The result Is 


Assuming the interstellar medium is mainly H, with a 10% admixture of He, 
we have <n?> = 1.3 amu. Taking the propagation energy to be E = 600 MeV/ 
nuc, and putting (10 - 2.2 x 10^ yr (Yiou et al. 1972), we have 


Using equation III. 42 for the total cross section, we also have 




(IV.19) 


10 . 


It now remains only to estimate the abundance of Be expected 
~ND 

if none decays, This can be done by inverting equation IV, 11 to 

calculate the expected observed abundances given an assumed set of source 
abundances : 


7J.r) = i \ Q^(y) 

^ 0 « 


(IV. 20) 


The source abundances used are, of course, based on the results of the 
previous calculation, when the leakage pathlength, \ , was estimated. 

This has been done by Tsao et al. (1973) at high energies (E > 2300 MeV/ 
nuc) using the. observed charge abundances and assuming the isotopic 

N 

composition of the various elements in the cosmic rays at their source is 
the same as that observed for solar system material, in. table IV. 3 we show 
the predictions of a similar calculation for the arriving cosmic rays, 
where we have adopted the source composition of Silberberg and Tsao, 



TABLE IV. 3: Compariforn With Model* 


Model 

Observed 

Z 

Y 

cm“^) 


5 

6 

Energy (lleV/nuc) 


600 

600 

Gaa Composition (percent He) 


10 

10 

Cross Section Modi£lcocions 


b 

b 

Abundance Ratios: 




’Be/Bo 

®Be/Be 

i“Be(Be 

0.61 -*0.07 
0.28 +0.05 
O.U +0.04** 

0.53 -*0.04 
0.28 +0.02 
0.19 +0.02** 

0.53 +0.04 
0.28 +0.02 
0,19 +3.02* 

Bc/C 

0.069+0.008 

0.090 

0.103 

1®B/B 

1»B/B 

0.30 +0.04 
0.70 +0.05 

0.32 -to. 03 
0.68 +0.06 

0.32 +0.03 
0.68 +0,06 

B/C 

0.251+0.017 

0.222 

0.248 

l»C/C 

ISC/C 

0.94 -*0.03 
0.06 +0.01 

0.95 +0.04 
0.05 +3.01 

0,94 +0,04 
0.06 +0.01 

C/C 

1.000+0.036 

1.000 

1.000 


0.43 -*0.09 
0.57 +0.07 

0.62 +0.03 
0.38 +0.05 

0.61 +0.03 
0.39 +0.05 

N/C 

0.271+0.024 

0.247 

0,260 


X 

U 

V 

u 


5 

5 

5 

5 


1200 

600 

600 

600 


10 

0 

10 

10 


b 

b 

a 

c 


0.54 +0.04 
0,26 +0.02 
0.20 +0,02** 

0.52 +0,04 
0.28 +0.02 
0.20 +3.02** 

0.58 +0.05 
0.23 +0.02 
0.19 +3.02** 

0.55 +0,04 
0.28 +0.02 
0.16 +0.01** 


0.101 

0.098 

0.055 

0.037 


0.31 +0.03 
0,69 +0.06 

0.32 +0.03 
0.68 ^.06 

0.30 -*0.03 
0.70 +3,06 

0.33 +0.03 
0,67 +0.06 

H 

%> 

0.218 

0.251 

0.222 

0.225 

0.95 +0.04 
O.OS ±0.01 

0.94 +0.04 
0.06 +3.01 

0,95 +0.04 
0.05 +0.01 

0.95’ +0.04 
0.05 +0.01 

1.000 

0,62 +0.03 
0,38 +0.05 

1.000 

0.61 +0.03 
0.39 +0.05 

1.000 

0.51 +0.03 
0.39 +3.05 

1.000 

0.62 +0.03 
0.33 +0,05 

0.243 

0.263 

0.243 

0.247 

3 ^ 


,* With Correction* ACDF (fee table III. 3) 

** -“Be Stable 

+ Cross Section ModltieoCloa Schemes: a - pure semlemplrical fonmila 

b - somlesipirlcal formila normalized to results of Llndstrom eC al. (1975) for C & 0 prliaarlea 
(see cable IIl.A) 

c - semicmpirical. foxinula notmolitcd to results of Lindatrom et al. (1975) for C £. O prlmaric* 

(ace table III. 4) and data given by Sllberberg et al. (1973) for lifl ’•“Be (correction 
factor 0.59), ’’B - ’■“B (2,33), and s®Fo -* -’■B (0.69) 

•H- Error Units on model calculations estimated assuming lOX uncertainties on oU total cross sections and 15X uncertainty on partial cross section* 
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and performed the same calculation at a lower energy (E = 600 MeV/nuc) 

which is more appropriate to the present application. Table IV. 3 shows 

the effects of varying various parameters involved in this calculation. 

_2 

The specific examples include a leakage pathlength ® g cm (rather 

than 5 g cm , an interstellar medium consisting of pure H (without the 

10% He admixture), a propagation energy E = 1200 MeV/nuc (rather than 

600 MeV/nuc), and the correction of the partial cross sections used for 

discrepancies with the experimental, values tabulated by Silberberg and 

Tsao (1973) in cases not covered by the data of Lindstrom et al. (1975a). 

Prom the talbe, it is clear that the relative abundances predicted for the 

Be isotopes are quite insensitive to the above variations. By comparing 

the predicted and observed abundances of ^^Be, we estimate that (55 + 21)% 

of the ^^Be survives in the cosmic rays. Although the error here is 

10 

estimated solely from the statistical errors on the observed Be abund- 
ance as assigned by the maximum likelihood fit, it would not be signifi- 
cantly increased if cross section uncertainties were included. One thus 
has 


5icC« 



-t-c2.0 



(IV. 21) 


If the leakage pathlength is = 5 g cm , we now have from equation 


IV. 17 


’>■ I”'' :«) 




-3 


(IV. 22) 


and, consequently. 



(IV. 23) 
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Since the leakage pathlength is a rather uncertain parameter (even within 
the present data, as may be seen from figures IV. 1 and IV. 2), the explicit 
dependences of n and on are shown in figures IV. 3 and IV. 4, respect- 
ively. As can be seen from these figures, the above conclusions are not 
substantially altered by reasonable fluctuations in the leakage path- 
length, 

C. Conclusions 

The results given by equations -IV»22 and IV. 23 lie between the 
extremes of the disk and halo models, as discussed in section A. This 
is to be expected. However, the result is substantially more consistent 
with the disk model than the halo. In fact, if 'there is a halo at all, 
the present result seems to indicate that it should be substantially 
flattened. 

This interpretation of the experimental data is obviously model 
dependent. In particular, in writing equation IV. 10 it has been assumed 
that the gas densities encountered by the cosmic rays along their trajec- 
tories through interstellar space is everywhere constant. If this is not 
the case, then the density of interstellar gas, n, should be replaced in 
the model by 


71 
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(IV. 24) 


Thus the density given by equation IV, 22 is not a simple average of the 
gas density over the confinement volume, but rather a weighted average 
with the cosmic-ray intensity as the weighting function. This weighted 
average reduces to a simple average in cases where there is no correlation 
between cosmic-ray intensity and gas density. 




Figure IV. 3. Matter density in cosmic-ray confinement volume as a 


function of mean escape pathlength. 


i 
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Figure IV. 4. Mean lifetime against escape for cosmic rays as a function of mean 
pathlength for escape. 


141 



-142 


Models which have been constructed to explain the gamma-ray emission 
of the galactic disk generally require a positive correlation between 
cosmic -ray intensity and gas density. Such a correlation would tend to 
make n somewhat larger than a simple average over the confinement volume 
of the gas density. Consider for example the model of Bignami et al. 

(1975), in which the cosmic-ray intensity is taken to be proportional to 
the gas density. In this model, the confinement volume is taken as the 
galactic disk, which is considered approximately equally divided between 
spiral arms and interarm regions. The gas densities in the spiral arms 
are taken as roughly a factor of two greater than those characteristic 
of the interarm regions. We then estimate that the weighted average gas 
density, IT, is only a factor of greater than the simple average of the 
gas density over the galactic disk. Thus it seems probable that n does 
not deviate significantly (at the current level of resolution) from a 
simple average of the gas density. 

If the relative composition of the cosmic rays (as opposed to inten- 
sity) is correlated with gas density, then the weighted average gas density, 
n, should also depend upon the cosmic-ray species from which it is deter- 
mined. This dependence would also be reflected in the pathlength = 

<nt> npcTg. Such variations in n, however, should be no more than those 
estimated above for intensity variations, and thus are almost certainly 
negligible for purposes of interpreting the present data. 

In general, the cosmic -ray source function used in the leaky box 
model, ^( 7 ), represents an average over space and time of the true source 
function, Q^( 7 ,r,t). Thus the true source function, which is almost 
certainly a highly discrete function of space and time, is replaced in the 
leaky box model by a constant in space and time. As we have already 



noted, fhis is generally a good approximation as long as none of the 
discrete sources are too close to the observer in space and/or in time. 
However, if a significant part of the pathlength, Xe» is traversed in a 
high density region near the source (e.g. a supernova envelope), we 
would expect to observe a smaller surviving fraction of ^®Be than if the 
density were constant everywhere along the typical cosmic-ray trajectory. 
Thus both the mean escape time, and the gas density in the correspond- 

ing confinement volume, n, would tend to be overestimated under such 
conditions . 

If the cosmic-ray electron component is assumed to be produced 
by the same sources and to be confined to the same "leaky box" as the 
nuclear component, then the mean escape time, Tg, should be reflected 
in a break observed in the electron energy spectrum. So far, no such 
bi'eak has been seen in observations extending from roughly 10 GeV (below 
which the spectrum is significantly affected hy solar modulation) to about 
ICOO GeV. The Interpretation of this observation is substantially clouded 
by the differences among the results reported by the various experimenters, 
which are significantly outside of the quoted experimental errors. 
0bser\<’ation3 of the electron -energy spectrum have been reported described 

as power laws with spectral indices ranging all the way from ~ 2,7 (which 

« 

would be consistent with the nuclear component) to ~ 3,4. Until the 
experimental situation is improved, it is not clear that observations of 
the electron spectrum can be used to place significant constraints on the 
mean confinement time of cosmic rays, Tq. Furthermore, such a result 
requires additional assumptions before it can be related to the parameter, 
Tg, which has been estimated in the present thesis. In particular, one 
generally assumes that the sources, confinement volume and mean escape 
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time of the electron component are identical to those of the nuclear 
component j none of which is established. For further details on this 
problem, see Meegan et al. (1975) and references therein. 

Let us now consider the resolution of the present experiment, 
and how it might be improved. The most obvious possibility would be to 
perform the same experiment at the top of the atmosphere. This would 
improve the uncertainty on the abundance at the top of the atmos- 

phere from 407o to ~ 25%. The improvement would not be as marked for other 
isotopes, however, since the fraction of the observed abundances of these 
isotopes originating in atmospheric interactions of heavier species is 
generally much smaller than for One might also hope to benefit 

from improved statistics due to the increased exposure time of an exper- 
iment done at the top of the atmosphere (i.e. a satellite experiment). 

The alternatives discussed above decrease the mathematical 
uncertainty in the result without a corresponding increase in the mass 
resolution of the experiment. Thus the credibility of the result may 
not in fact be significantly improved. The most important possibility 
for improving the actual mass resolution of the experiment would be an 
improved calibration of the various detectors. Such a calibration would , 
most probably be performed at an accelerator and include the internor- 
malization of the detector responses, a detailed investigation of the 
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shapes of the response functions of the various detectors for various 
particles of various incident velocities, and the mapping o'f positional 
variations in detector response. A knowledge of the shape of the response 
function is particularly important if we are to improve the credibility 
of the present result; specifically, the question of a possible asym- 
metry in the mass distributions for Be isotopes urgently requires 
resolution. The mass resolution may also be improved by as much as a 
factor of two by going to completely unsaturating scintillators (if such 
can be found). This would' also help to improve the situation with 
respect to the internormalization of the detectors, since the response 
functions of such scintillators would be expected to be more nearly 
uniform. We see from the above discussion that there is a large potential 
for achieving improved mass resolution from experiments of this type. 

Furthermore, only a small part of this potential need actually be 

9 10 

achieved before a valley begins to separate the Be and Be distri- 
butions, as discussed in Chapter II. Such a situation would, of course, 
markedly improve the credibility of the present results. 

Ideally, future experiments should be done at higher energies in 
order to minimize the importance of corrections for solar modulation. 
Energies ^ 1 GeV/nuc would be most advantageous if achievable. One 
possibility for an isotope experiment operating at such energies would 
be to measure the parameters rigidity, Cerenkov radiation and as 
was suggested in Chapter II. This has the advantage of not requiring 
that all a particle's energy be absorbed before the measurement can be 
carried out. The correction for interactions in the detector is accord- 
ingly less severe. On the other hand, experiments where the partiele 
must be stopped in the experiment suffer increasing losses to interactions 
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in the detector as particle energies increase, and their statistical 
weight is accordingly reduced. In addition, all classes of experiment 
must contend with energy spectra which are rapidly falling with increas- 
ing energy. Thus, experiments to directly measure the abundances of 
rare isotopes such as ^*^Be in the cosmic rays at high energies (E > 

1 GeV/nuc) may not be practical without long exposure times such as 
might be obtained in satellite experiments, for example. 



APPENDIX A 


DETECTOR CALIBRATION 

This appendix deals with the detailed calibration of the instru- 
mental response using in flight data. In order to do this, it is re- 
quired to determine the gain and zero offset and associated 

with each gain range of each pulse height analyzer (i.a. see’ equation 
III. 8), and the correction function (M^) for positional and temporal 
variations in detector response. 

1. Gain Change Factors . The output of the pulse height analyzers 
is an integer pulse height together with a two bit code to indicate the 
appropriate gain range. The use of automatic gain switching was re- 

5 

quired in order to obtain a dynamic range of 10 . This necessitated 
the reconversion of the resulting pulse heights to a linear scale before 
the reduction of the data could proceed further. Assuming the pulse 
height analyzers and amplifiers are linear in all gain ranges, we recall 
equation III. 8: 


H. = ( N j -t 


(A.l) 


Here is the final, linearized pulse height characterizing the response 
of detector i, and is the same quantity as defined by equation III. 2. 

denotes the unnormalized pulse height from the gain range'of 
detector i, while and b^j refer to the gain and zero offset required 
to normalize it, respectively. The usual procedure would be to put the 
b. . = 0, then observe that at a gain change boundary the same pulse 

J 

height is measured twice, once in each of two different gain ranges: 




(A, 2) 
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TABLE A.l: Gain Chang« Factors 





DETECTOR 

1 

2 

3 

4 

5 

6 

7 

1 

2 

3 

4 

5 

6 7 


D1 

0.142 

3.96 

4.14 

4.00 

— 

— 


2.99 

0 

0 

0 

— 



D2 

0.179 

8.58 


H m sins 

M n M «• 

4»«**S* 

as • 

0 

0 


— 

— 



ll. 

D3 

0.142 

0,389 

3.13 

23.4 

0.188 

0.755 

.3.02 

3.30 

-3.86 

1.61 

3.03 

-32,82 

-4,20 

C* 

D4 

0.394 

2.55 

2.40 

2.52 

2.16 

2.50 

— 

8.48 

2.14 

-1.99 

1.07 

1.89 

-0.76 

D5 

0.420 

2.24 

2.72 

2.22 

2.39 

2.25 

— 

10.05 

1.89 

-1.67 

0.64 

1.89 

-0.84 

D6 

0.367 

2.34 

2.56 

2.28 

2.25 

2.25 

— 

8.31 

4.63 

-1.75 

1.47 

1.78 

-0.79 

D7 

0.0923 

2.37 

2.44 

2.42 

2.09 

2.35 

— 

16.14 

1.26 

0.27 

2.81 

2.89 

-1.23 


D8 

0.0983 

2.30 

2.63 

2.28 

2.23 

2.31 


19.13 

1.84 

-1.15 

0.67 

3.37 

-1.46 


D9 

0.1117 

2.45 

2.43 

2.39 

2.05 

2.11 

f" 

3.48 

3.68 

0.23 

0.13 

5.12 

-2.43 


DIO 

0,1040 

2.28 

2.68 

2.20 

2,34 

2.16 

— 

13.27 

1.65 

-1.23 

0.56 

2,85 

-1.32 


DIL 

0.1037 

2.30 

2.53 

2.47 

2.09 

2.38 


10.13 

0.17 

1.28 

-1.36 

2.57 

-1.08 


D12 

0.1070 

2.34 

2.67 

2.21 

2.30 

2.09 


17.66 

2.86 

-2.00 

1,02 

2.66 

-1.27 — - 


D13 

0.0986 

2.48 

2.42 

2.52 

1.86 

2.60 


7.30 

2.06 

-0.79 

0.54 

5.00 

-1.92 


D14 

uO.0942 

2.18 

2,89 

2.04 

2.58 

1.99 

.... 

6.40 

2.16 

-0.86 

0.65 

2.42 

-1.22 



JJ- 
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i>\ 


(A. 3) 


In the above, we note that gain ranges are numbered with decreasing gain. 
The a., for detector i are thus all determined in terms of a.,, which is 

ij - 2.1 » 

recognized as the normalization constant for detector i, -and will be 
determined elsewhere. This is the procedure actually used for detectors 
Dl, D2 and for the highest gain change of D3 (where the 256 channel 
analyzer has saturated, and no longer overlaps the 1024 channel analyzer). 
For the 1024 analyzers, the a^^ are thus determined to a maximum pre- 
cision of about 0.5%. If the zero offsets, b--, are of the order of 5 

J 

channels , then additional errors of the order of ~ 2% could be intro- 
duced at each gain change. The a,, resulting from the analysis just 
described, as applied to detectors Dl, D2, and the last gain change of 
detector D3 are tabulated in table A.l. 

A substantial innovation included in the current experiment was 
the use of dual pulse height analyzers with staggered gain ranges on 
detectors D4-14 and, to some extent, on D3, This offered the advantage 
of being able to determine both the and the much more accurately. 
In general, than would be possible using the method just described. In 
addition, this made it possible to detect and correct for noise, bit 
errors and other problems related to the amplifiers and pulse height 
analyzers at a very early stage of the data reduction. If we again 
order the gain ranges (from either pulse height analyzer) with j increas- 
ing as- gain decreases we have, from equation (A.l), 



r 


U'' 







(A.4) 
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(A.5) 


for the region where the gain ranges j and j+1 overlap. Thus the gains, 

a , and zero intercepts, h- can all be determined in terms of a-i 
ij J 1 

and which are recognized as the absolute normalization and zero 

intercept parameters for detector i (as a whole), and will be determined 
later. Figure A.l shows a representative plot of the type just described, 
in this case for and most cases, we feel the relative gain 

factors, a^£j,have been determined to within — 0.1%, while the zero off- 
sets, b.., have been determined to within ~ 0,2 channels by this approach. 

J 

The detectors D8 and D9 show the presence of substantial amounts of 
noise, probably originating in the spark chamber, and are accordingly not 
as well calibrated. Figure A. 2 shows versus which is by far 

the worst case. In addition to noise, diagrams such as figure A.l are 
sensitive to bit errors, which appear as points well off of the cali- 
bration line. Nohlinearities' in the amplifiers can also be detected 
and corrected. Finally, the use of an additional pulse height analyzer 
for each detector increased the resolution immediately after a gain 
change by a factor of about 2,3 over the case where only one pulse 
height analyzer is used for each detector. This follows since the 
resolution is determined by the largest of the two pulse heights which 
(after the first gain change) is never less than about 90; for the same 


signal, the smaller pulse height (i.e. the one which has just gain changed) 
is about 40. 

The pulse heights from the two pulse height analyzers for detector 
i, nr. and were combined to obtain a single value, H., according 


to the relation 
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Figure A.l. 


Typical gain change calibration plot: 
versus 
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0-4 -f a,^-,. 

Here it is assumed the resolutions of the pulse height analyzers (in 
channels) are given by 


so that we have 


() f Consfcun't 

H 


(A. 7) 


S = 


(A. 8) 


from equation A.l. This is equivalent to taking the normal from the 


experimental point ^i('+l)^ appropriate calibration line, 


The resolution of is given by 


- >V [^i + ^44 " S o.U 07;, 


(A.9) 


where a Sf 2.4 has been assumed. Putting o ^ ^ 0.5 channels, we 

i(j+l) H 

have 


0- < 0.46 0 .^ 




IPW.X 


(A. 10) 


where j is the gain range giving the largest pulse height, HT., 
max xj 

corresponding to 

Events which lie clearly off th« calibration are tagged and, 
usually, thrown out in the analysis. This is intended mainly to remove 
events with bit errors. The resolution for doing this is generally a 
factor of about 2.4 worse for the larger of the two pulse heights, 
because of the slope of the calibration line. 

For detector D3, all but the last gain change of the 1024 channel 
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analyzer were calibrated in the fashion just described* Since the last 
gain range of the 256 channel analyzer had saturated before the last 
gain change of the 1024 channel analyzer was reached, this last gain 
change factor was determined by the same procedure as for Dl and D2, 

The values of a^- . and . obtained by the above procedure for detectors 

* J J 

D3-14 are tabulated in table A,l, 


2, Detector Nonlinearities . The amplifiers used (especially-' 
with the 256 channel analyzers) were known exhibit some nonlinearities, 

A pulser calibration of a representative amplifier - pulse height 
analyzer combination is shown in figure A. 3. It is‘ apparent from the 
figure that, above channel 32 the response is nearly linear with the 
nonlinearity being manifested mainly in a large zero offset. For this 
reason, the calibration lines (e.g. figure A.l) were fitted neglecting data 
below channel 32 in either pulse height analyzer. The nonlinearities 

have otherwise been neglected. This can only significantly effect the 
analysis for extremely low pulse heights in the first gain range (i.e. 

< 15 X minimum for D4-6 and < 5 x minimum for D7-14). 

3. Detector Normalization . Originally it had been intended to 
use highly relativistic carbon and oxygen triggering the experiment in 
the penetrating mode to determine the detector normalization factors. 


a and to provide a first order estimate of the zero offsets, b.,, 
il il 

This could be done, for example by selecting = 0, then making a plot 
of the relativistic peaks (times minimum units) of various elements in 


Dl against their expected (after normalization) values for an ideal. 


non- saturating scintillator. Because of scintillator saturation (non- 
linearity) effects which increase with increasing charge, the result would 
be a curve asymptotic to a straight line through zero in the limit Z -> 0. 



CHANNEL NUMBER OUT 


260 


REPRESENTATIVE 256 CHANNEL 
PULSE HEIGHT ANALYZER 
LINEARITY CALIBRATION 



SIGNAL AMPLITUDE 


Figure A. 3. Representative calibration of nonlinearity 
of a 256 channel analyzer with a pulser. 
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The slope of this line would then be identified with the parameter 
The remainder of the scintillation detectors, D3-14, would then be 
easily normalized by adjusting the and b.. so that both relativistic 

carbon and oxygen fall in the same channel as for Dl. The Cerenkov 
detector was indeed finally normalized on this basis (i.e. by setting 
relativistic carbon equal to 36 x minimum and choosing ~ > however, 

it was found that the above procedure could not be directly applied to 
the normalization of the scintillation detectors. 


When the various scintillation detector pulse heights (unnormal- 
ized) were plotted against the Cerenkov, it was noticed that the scintil- 
lator response to the most highly relativistic particles increased as 
one went deeper into the stack. This was attributed to very high energy 
delta rays which penetrated several detectors. Thus it was decided to 
use for the normalizations semi-relativistic particles whose velocities 
were not sufficient to give, rise to such high energy delta rays. In 
order to use these particles, however, one must correct for the slowing 
down of the particles as they penetrate into the stack. This was done 
by first determining a particle's range from the top of the Cerenkov 


detector, R 2 J from equation III. 6; 

^ C dR' 



(A. 11) 


Here the range -energy tables (Barkas and Berger, 1964) are used to 
relate the range to the velocity, together with the assumption that mass 
and charge , are related according to A =* 2Z (this is excellent for 'C, 0, 
"Mg and Si; in any case, the corrections are small enough that the 
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departures from this relation enter as second order terms, and may be 
neglected). The approximation in equation A. 11 is that the particle 
does not slow down significantly in traversing the Cerenkov radiator, D2; 
this is not required for the analysis. The parameter C' is unity for 
highly relativistic particles, and zero for particles which emit no 
Cerenkov light. If the rate of radiation of Cerenkov light changes 
significantly over the Cerenkov detector, the solution of A. 11 is by an 
iterative procedure. 

The quantity 





Sec6 






(A. 12) 


is now evaluated for each detector, i. The particles are selected for 
the analysis according to 0.1 ^ CT ^ 0,8, which is essentially the 
restriction to semi-relativistic particles mentioned above. It is then 
noted that a plot of against 1/2 (Hj^ + H^) shows a very nearly linear 
behavior (for a given charge) in this region, as is shown in figure A. 4: 



P 

Z 


( ^j) 



(A. 13) 


It is assumed that and C(- should be related by exactly the same co- 
efficients ( )) and ) if the detectors are all properly normalized, 

i 1 e • 


cl = PH.tt 


(A. 14) 


Using equation A. 14, the a,^ and b^j^ for all scintillation detectors 
may be related to those of one standard detector by plotting 



(A. 15) 



I/2(H, + H3)(xM1N) 
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Hg(xMiN) — ^ 

Figure A. 4. Unmapped Cerenkov versus scintillator plot. Pulse 
heights are in times minimum units such that a 
particle of velocity g=l and charge Z=1 gives a 
pulse height of 1, regardless of detector thickness 
or angle of incidence. 
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against standard detector gives the optimal values for the 

parameters j) and t . As the analysis is described above, this role would 
be played by the combination of .Dl and D3 (see equation A. 13). However, 
because of the unique response functions characterizing these two 
detectors, it was preferred to use D4 as the standard detector. The 
normalization constant was determined from the relation (see 

equation III. 3). 


d.L 

dz 




Cos 9 


Tr 






(A. 16) 


(A. 17) 


where the average is over all data with ^ 0.8 and 

h 0.8 

has been assumed for the normalization of the mapping function, The 

normalization is then obtained from the condition 


so that 


2^0 


db 

dE 


= I 


(A. 18) 


a 
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t: 

2-^0 




(A. 19) 




if we choose b =0. This is shown in figure A, 5. As is evident from 
41 

the figure, the missing data from the light (1<Z^3) nuclei (which should 
have been admitted in the calibration mode) would have been very useful 
for this purpose . 

Because of scintillar non-linearities which varied from detector 


•to detector, it was found that the interdetector normalization plots 
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just described did not always display linear correlations over the 
entire range of charges. In fact, the deviations from linearity became 
marked for highly charged particles, and especially for iron. It was 
therefore decided to use only particles with Z < 10 in the normalization. 

The normalization of detector responses obtained by the ■ procedure 
just described are throught to be accurate to within 2% for most cases. 
This is determined primarily by the available statistics; consequently, 
the accuracy is reduced for detectors deep in the stack where there are 
less particles available for the calibration, 

4, Positional and Temporal Nonuniformities . The variations in 
response of each detector (except the Cerenkov) as a function of position 
were corrected using a polynomia.l fit to the response of the detector 

to relativistic (C-' s 0.8) carbon and oxygen nuclei. If we now assume 

^ T 

that is separable into a time dependent part, M^(t), and a 

position dependent part, we have from equation III. 3 


< CJ>of 


(A. 20) 


One then obtains 


li^ = o.s 


(A. 21) 


T. 


cj > 0.8 


where 




(A. 22) 


has been assumed. It has also been assumed that the relative fluctuations 


in M, and T. are small. 
1 1 


P 

The function is approximated by an 11 x 11 polynomial in x^ 
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and y . 




A 

e. 


and y£ = 





9 


Mf = ^ ? :!■ £ yfe 


(A. 23) 


whose coefficients are determined by a tnaximuni likelihood fit. 

Since there are few experimental points near the edges of the scintil- 
lators, an imaginary border was introduced around each scintillator with 

p 

dumn^ data in order to prevent the rapid divergence of near the edges. 

The number of points available for fitting (equation A. 23) decreases by a 

factor of 2.5 as one goes from D1 to Dl4. Hence the accuracy of the map 

obtained for a detector decreases with increasing depth of the detector 

p 

in the stack. The resolution of the fit (equation A. 23) to thus 
varies from cr ^ 0.4% for Dl to cr =- i.5% for D14, Figures A. 6 and A. 7 


present conteur diagrams of the response of Dl to relativistic carbon and 
oxygen nuclei (x minimum units) , and the map fit to the data, respectively. 
The effectiveness of the map in correcting the Dl response may be seen 
by a comparison of figures A. 4 and A. 8. 

Slow drifts of the responses of the various detectors were dis- 
covered by observing the relativistic carbon and oxygen peaks as functions 
of time. It was found that the fractional deviations of oxygen and carbon 
peaks from their optimum values- were the same, so that there was no need 
to distinguish the two cases. This meant that there were no (observable) 

zero offsets involved, and we had only to determine gain variations of 

T 

the various detectors with time. This is corrected by the factor 
which is normalized so that 


<mT> = 1 


(A. 24) 









Figure A»/. Contour plot of Dl map, fitted to data of figure A. 6 
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in order to preserve detector normalizations. The function M, is then 

1 

determined by 




(A. 25) 


T 

The are linearly correlated, as is shown in figure A. 9 for the case 
of M^(t) andH^Ct). One thus has 

mJ *: o(J h 7 /3j ■ ■ (A.26) 

T T T 

and all the H.(t) are determined in terms of M, (t), the coefficients i?i • 

1 1 3 . 

T 

and being determined by least squares fits to data such as that of 
figure A. 9, Using theo(^ and 6^, is then approximated by a least 
squares spline technique (Thompson 1973) using equation A. 25 and data 
from all detectors. The result is shown in figure A. 10, together with 
the corresponding temperature profile in the gondola. There seems to be 
some correlation between the two. The normalization condition' (equation 
A. 24) is apparently not satisfied in this figure since the initial phases 
of the data reduction were carried out using only a fraction of the 
available data which was taken near the end of the flight. The normal- 
izations were then all based on this sample of the data, even when the 

analysis was expanded to include the entire data set. Within experimental 

T 

error, the function was parameterized using a least squares spline 
fitting routine (Thompson 1973). Since fluctuations in the data- were in 
general very much larger than the effects we were trying to measure, the 
data supplied to the fitting routine consisted of ~ 10 miniute averages 
of the actual data. 

Because of the statistics involved, we have been able to 
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[l- M^Ct) ] 100 


T T 

Figure A. 9. Linear correlation of drift functions and . 
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accurately determine only on time scales greater than about one 
hour. Although it is conceivable that variations in response on shorter 
time scales may be present (such as might be expected from a rotation 
of the gondola in the earth's magnetic field, for instance), we have no 
way to see or correct for these. 

5. Detector Thicknesses . The thickness of each detector was 
measured on a two inch grid after the flight. The thickness function, 
was then approximated by linear interpolations between these points. 
This correction was necessary because of the large fluctuations in the 
thickness of certain detectors--up to 10% in extreme cases. As an example, 
we have listed in table A. 2 the thicknesses of the various detectors at 
each of their four corners. We feel the detector thicknesses are known 
to ~ 0.1%. 
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TABLE A. 2 


Thicknesses of Detectors D4-14 at Coimers (in.) 



A 

B 

C 

D 

D4 

0.388 

0.352 

0,356 

0.377 

D5 

0.377 

0.361 

0.355 

0.395 

D6 

0.366 

0.361 

0.355 

0.371 

D7 

1.344 

1.431 

1.367 

1.338 

D8 

1.442 

1.362 

1.366 

1.388 

D9 

1.429 

1.345 

1.353 

1.351 

DIO 

1.356 

1,434 

1.381 

1.371 

Dll 

1.284 

1.335 

1.348 

1.351 

D12 

1.325 

1.376 

1.369 

1.498 

D13 

1.387 

1.390 

1.399 

1.615 

D14 

1.388 

1.619 

1.395 

1.374 



APPENDIX B 


aij 




DEFINITIONS OF SYMBOLS USED IN THIS THESIS 

= gain factor for gain, range of detector i, 

= gain factor for gain range of detector i relative to 

(j “ 1)^^ gain range. 


A t= mass number. 


Ai = mass deduced from the experimental point (L^jRi) . 

Ai(6.) = mass as a function of velocity and charge as deduced from the 

observation of one of the observables i = E, R, P or C. 

Ag = mass characterizing a standard response curve. 

b^j = zero offset for gain range of detector i. 

— ij ~ zero offset for gain range of detector i relative to 

(j “ 1)*^^ gain range. 


B = magnetic field intensity. 

gauss = magnetic field intensity characteristic of the 
galactic disk. 


Bj^ 10 gauss = magnetic field intensity characteristic of the 
galactic halo. 

c = 2.998 X 10^® cm sec ^ — speed of light. 

C* ~ Cerenkov pulse height. 




A-Atmos 


A^E 

Z 


A^nlnter 

Z 


*D 

LBj 

^D 


value of expected if detector i is actually .a Cerenkov 
detector of the same thickness, 

atmospheric correction factor for particles of mass A and 
charge Z. , • 

correction factor for energy windows for particles of mass A 
and charge Z. 

correction factor for losses to interactions in the detector 
(including the gondola) for particles of mass A and charge Z. 

atmospheric propagator matrix (= weighted slab propagator 
matrix) . 

leaky box propagator matrix, 
slab propagator matrix. 
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§ = unit vector orthogonal to directed towards electronics 

package, 

= unit vector pointing towards zenith. 

= ( 7 -I) mc^ = kinetic energy of a particle of mass A, charge Z 
^ and velocity g. 


= ionization energy loss or stopping power of a particle of 
^ ^ mass-Aj charge Z and velocity p. 


£(X) = distribution of atmospheric depths. 


= mean free path for production of species j from spallation 
of species i. 


= geometrical factor. 


= mean free path for destruction of species i. 




calibration pulse height for detector i as predicted from 
Cerenkov and d4. 


= energy deposited in detector i by particle of mass A, 

■■ charge Z and velocity p, 

if. = (integer) pulse height from gain range of detector i. 


= ratio of light observed in last two detectors to that 
observed in last detector. 

12. 5Z eV = average ionization potential of atom of atomic 
number Z . 

= differential energy intensity of species i. 


over space. 


= differential rigidity intensity of species i. 


= flux of species i observed under a slab of thickness X. 
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J,.(X) = 
ik 


flux of Species i observed under a slab of thickness X 
which has undergone exactly k interactions. 


= diffusion coefficient. 


28 


s: 10 


Kjj ^ 1Q29 


Kp 

K 


9 —1 

cm^ sec"'*' = diffusion coefficient characteristic of the 
galactic disk. 

9 —1 

cm sec = diffusion coefficient characteristic of the 
galactic halo. 

rigidity dependent part of diffusion coefficient K=8Kj.Kj,. 
radial dependent part of diffusion coefficient K=pKj.Kp. 


\(P) = 

<7 i 




total light observed due to particle of mass A, charge Z and 
velocity p coining to rest in detector i. 


= differential saturation function for response of scintillator 
i to a particle of charge Z, mass A and velocity 8* Usually 
taken to be independent of mass and detector number. 


m 


= mass. 


n = 9.11 
e 


OO 

X 10 g = electron mass. 


= 1.672 X 10~^^g = proton mass. 


X minimum unit = pulse height expected from a relativistic (p=l) particle 
of unit charge and mass which traverses the detector 
along the same trajectory as the particle in question. 

= factor describing positional and temporal variations in 
detector response. 


1 X 


P T 

= position dependent part of Mi. 


T 

M^(t) 


= time dependent part of M. = mT 


1 1 


n 


n 


n Si 
D 


= number density, 

= index of refraction. 

“3 

1 atom cm = atomic number density characteristic of the 
galactic disk. 
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= density of cosmic rays of species i 


-2 ”3 

Tijj ai 10 atom cm 


= atomic number density characteristic of the 
galactic halo. 


’N. 


tot 


= total number of events observed corresponding to mass A and 
charge Z. 


N(A"*) = number of events observed in dA about A corresponding to 

Z an isotope of actual mass A- and charge Z. 


2i 


= (numerical) observed value of the observables i = E, 
P or C. . , 



,0i(p) 


= functional dependence of one of the observables i = E, 
R, P or C on Z, A and g. 


dE 

dx’ 


p = yprac = momentum. 

” power of the cosmic-ray sources. 

QR.S 



= rigidity of a particle of mass A, charge Z and velocity 0. 



r 


= source function for species i (including nuclear interactions) 

= source function for species i (excluding nuclear interactions) 
s 

= averaged over space. 

= heliocentric radius. 


r 

c 



= gyroradius of a particle of momentum, p, and charge Z, 
in a magnetic field, B. 


-13 

r = 2.818 X 10 cm = classical electron radius, 
e 

r„ S: 12 kpc = radius of galactic halo. 

H 

= positional vector with respect to the geometrical center of 
detector i. 
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= range tneasured from the top of detector i. 

= range of a particle of mass A, charge Z and velocity g. 

= surface vector (orthogonal to surface, pointing out). 

= time. 

Tj^(r^) = thickness of detector i as a function of position. 

V 21 400 km sec ^ - solar wind velocity. 

V = volume 

-3 

w 21 1 eV cm = local cosmic-ray energy density. 

21 250 pc = half thickness of the galactic disk. 

Z = charge or atomic number. 

a = coefficient 

T T T 

= linear correlation coefficient (slope) for versus 



p = velocity relative to that of light, 

X _ T ■ T 

= zero intercept of linear approximation to versus 


Po’^oi “ intercept parameter in linear approximation to A^(p) . 
2 - 1/2 

7 = (1 - p'^) = relativistic Lorentz factor. 


. = mass resolution scale factor for contribution to total ' 

resolution from resolution of the measurement of the 
observable i = E, R, P or C in a two parameter 
analysis . 

6 = exponent in power law approximation to range-energy relation, 

10 ^ = cosmic -ray anisotropy. 

CR 

r| = zero intercept of versus § (H^^ + H^) . 

0 = zenith angle (with respect to detector normal). 
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-2 

3: 5 g cm = mean pathlength traversed by the cosmic rays. 


= characteristic dimension of a luminescence center, 


V = slope of versus + H 2 ) . 

f — mass density. 

O' = resolution in terms of the standard deviation characterizing 

an (assumed) gaussian distribution. 


= resolution of a pulse height analyzer (in channels). 


'H 


= energy resolution of the gain range of detector i due 


ij to pulse height analyzer round off errors, 




partial cross section for production of species j from 
interaction of species i with species k. 


= total inelastic (reaction) cross section for 
interaction of particles of masses A and 


T^(j -»■ i) = mean lifetime against radioactive decay taking species 
j to species i (* denotes "anything")* 


Cf = azimuthal angle (about zenith). 

= streaming of cosmic rays of species i. 


yCR 
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